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FOREWORD 


The space shuttle main engine (SSME), a reusable space propulsion system. 
Is the most sophisticated Earth-to-orblt propulsion system In use today. Since 
the United States will rely on the shuttle and Its derivative versions for 
space transportation for tfie next several decades, the National Aeronautics and 
Space Administration (NASA) Is conducting a program which will establish a 
technology base and develop the analytical tools necessary for an orderly evo- 
lution of reusable rocket engines. 

A substantial portion of the technology program Is devoted to the struc- 
tural Integrity and durability aspects of reusable rocket systems. This effort 
focuses on the development of accurate analytical models to describe flow 
fields, aerothermodynamlc loads, structural responses, and fatigue/fracture 
effects. These analytical tools are then applied to the development of life 
prediction models of engine components and subsystems. The effort Is supported 
by the development of advanced Instrumentation with capabilities to verify the 
models In an SSME-Ilke environment. These Instruments have the potential for 
future use as diagnostic sensors to perform real-time condition monitoring In 
critical engine components. The research Is conducted through In-house activi- 
ties at the Lewis Research Center, Industrial contracts, and university grants. 

To provide a forum for representatives from government. Industry, and 
academia to learn about and discuss the latest findings and progress toward 
Improved structural Integrity and durability of reusable space propulsion 
systems, a 2-day conference was held at NASA Lewis In May 1987. This publica- 
tion contains extended abstracts and figures used for each of the presentations 
at the conference. 


Stanley 0. Marslk 
Conference Chairman 
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ADVANCED EARTH-TO-ORBIT PROPULSION TECHNOLOGY 


Frank W. Stephenson 
NASA Headquarters 
Washington, D.C. 


The Structural Integrity and Durability Program Is an Integral part of a 
broadly based NASA rocket technology effort addressing reusable Earth- to-orblt 
(ETO) engine Issues which Is being conducted through the NASA Marshall Space 
Flight and Lewis Research Centers. Program funding Is provided by the Office 
of Aeronautics and Space Technology (OAST) and the Office of Space Flight 
( OSF ) . 

The era of reusable rocket engines was ushered In by the space shuttle 
main engine (SSME). A small effort aimed at extending reusable engine technol- 
ogy was Initiated In the 1970 * s , but It wasn't until 1980 that full attention 
was given to this critical area. At that time, it was fully recognized that 
our ability to characterize, predict, and guarantee the life of reusable rocket 
engines, such as the SSME, was severely limited. Since that time, a focused 
program aimed at developing design and analysis tools capable of assuring the 
life and performance of reusable ETO engines has been strongly supported on a 
continuous basis. The program has used SSME experience to Identify critical 
life-limited engine components, as well as design and analysis deficiencies, 
and In part Is building on the experience base established In the airbreathing 
engine community to help correct these deficiencies. The SSME will also be 
used as a test bed to verify advanced design concepts and to validate Improved 
analytical codes developed under this program. 

The ETO program encompasses three major areas: life prediction/extension, 

performance, and operations. The life prediction/extension efforts are focused 
on developing an understanding of and analytically simulating Internal engine 
environments, both mechanical and thermal; on defining resulting transient and 
steady-state loads; on developing analyses for determining component structural 
response and material behavior; and on ultimately developing life-prediction 
techniques, especially for combinations of high-cycle and low-cycle fatigue. 
Engine performance technologies are primarily directed toward extending combus- 
tion, stability, heat-transfer, and cooling models to higher pressure operating 
conditions. Higher temperature turbine blade materials are also being explored 
to allow higher combustion pressures. Operations technology Is directed toward 
developing health monitoring diagnostic sensors for automated between-f 1 Ight 
Inspection, servicing, and checkout; for maintenance scheduling; and for In- 
flight fault-tolerant engine operation. The principal beneficiary of these 
technology advances will be the next generation of ETO rocket engines, although 
the SSME has already accrued some benefits. This conference Is focused primar- 
ily on reviewing the status and results of life extension/prediction technol- 
ogy tasks Involving engine hot-gas flow systems. 
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OVERVIEW OF AEROTHERMODYNAMIC LOADS DEFINITION STUDY 


Raymond E. Gaugler 
NASA Lewis Research Center 
Cleveland, Ohio 


The objective of the Aerothermodynamlc Loads Definition Study Is to 
develop methods to more accurately predict the operating environment In the 
space shuttle main engine (SSME) powerhead. Development of time-averaged and 
time-dependent, three-dimensional viscous computer codes as well as experi- 
mental verification and engine diagnostic testing are considered to be essen- 
tial In achieving that objective. Time-averaged, nonsteady, and transient 
operating loads must all be well defined In order to accurately predict power- 
head life. Improvements In the structural durability of the SSME turbine drive 
systems will depend on our knowledge of the aerothermodynamlc behavior of the 
flow through the preburner, turbine, turnaround duct, gas manifold, and 
Injector post regions. 

The approach taken In this study consists of two parts: (1) to modify/ 

apply/dlssemlnate existing computational fluid dynamics (CFD) tools In response 
to current needs, and (2) to develop new technology that will enable more accu- 
rate computation of the time-average and unsteady aerothermodynamlc loads In 
the SSME powerhead. With the more accurate aerothermodynamlc loads predictions 
providing boundary conditions to Improve structural and fatigue life analysis, 
the goal of Improved durability will be met. 

This study was begun In October 1983. The Initial effort Involved the use 
of existing CFD tools to analyze problems In the fuel and oxidizer turbopump 
turbines, the fuel turbine turnaround duct, the fuel-side preburner, and the 
main Injector liquid oxygen posts. Results of those efforts have been pre- 
sented at previous meetings. 

At this meeting, some of the currently ongoing tasks will be highlighted, 
primarily those dealing with the development of new CFD tools. Also Included 
will be a final report on the analytical and experimental study of the main 
Injector liquid oxygen posts. 
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UNSTEADY STATOR/ROTOR INTERACTION 

Philip C.E. Jorgenson and Rodrlck V. Chlma 
NASA Lewis Research Center 
Cleveland, Ohio 


The major thrust of the computational analysis of turbomachinery to date 
has been the steady-state solution of Isolated blades using mass-averaged Inlet 
and exit conditions. Unsteady flows differ from the steady solution due to 
Interaction of pressure waves and wakes between blade rows. To predict the 
actual complex flow conditions one must look at the time accurate solution of 
the entire turbomachine. 


ANALYSIS 

Many of the numerical tools used In the analysis of Isolated blades can 
be used for time-accurate analysis. Here the quasl-three-dlmenslonal Euler and 
thin- layer Navler-Stokes equations are solved for unsteady turbomachinery 
flows. These equations are written In general coordinates for an axlsymmetrlc 
stream surface, and they account for the effects of blade- row rotation, radius 
change, and stream- surface thickness (ref. 1). 

A four-stage Runge-Kutta scheme based on the work of Jameson (ref. 2) Is 
used to predict time-accurate results. Body-fitted C-type grids were used In 
this work and were generated using the GRAPE code (GRIds about Airfoils using 
Poisson's Equation) developed by Sorenson (ref. 3). 

To predict the Interaction between a stator and a rotor In turbomachinery, 
one must pass flow Information to an Interface that acts as a moving boundary 
between the two computational grids. A nonconservative Interface formulation 
Is used In this procedure. The solution must be Integrated using a constant 
minimum time step based on the computational domain of two blades. The data 
management necessary to update the stator and rotor flowflelds will be dis- 
cussed. Currently the code Is limited to solving stator/rotor configurations 
that have equal pitch blades. 


RESULT S 

The solution procedure has been applied to two test problems. The first 
Is a cascade and the second Is a turbine rotor from the space shuttle main 
engine (SSME). Euler and Navler-Stokes results will be presented. 

A model turbine stage made up of two Identical cascades of NACA 0012 air- 
foils was used to develop and test the Interface and data management routines. 
Figure 1 shows the C-grlds used In the computation where the second blade row 
Is moving downward. A converged solution can be seen In the periodic loading 
diagram of the rotor (fig. 2). 
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The stator Inlet Mach number Is 1.84. Figure 3 shows relative Mach number 
contours when the blade rows are aligned. Detached bow shocks form ahead of 
the stators. Interact at mid-passage, and reflect obliquely back to the stator 
surface. A strong curved shock forms at the stator trailing edge, and the 
shock curvature generates an entropy layer that convects downstream into the 
rotor. The stator flowfield Is symmetric but the rotor flowfleld is asym 
metric due to the incidence of the relative flow. 

Figure 4 shows relative Mach number contours after the rotor has moved 1/4 
pitch downward. The flow between the stator and rotor is subsonic so that 
pressure waves from the rotor can affect the stator trai 1 ing- edge shock, which 
has moved upstream. The entropy layer from this shock produces a large asym 
metry In the rotor flowfield. lhe flow reaccelerates to supersonic in the 
rotor passage and produces strong oblique shocks off the rotor trailing edge. 

Previous steady-state calculations of the SSMF first-stage fuel- side tur- 
bine blade (ref. 1) showed a reverse flow region on the pressure surface 
(fig. 5). lime- accurate results are being generated to determine how that 
separation is affected by the upstream stator wakes, lhe grids used for this 
calculation are shown in figure 6 where the second blade row is the rotor of 
Interest. Euler results are shown in figure 7. There Is very little interac- 
tion between the blade rows. Viscous results will be shown in the presenta 
tion. Currently the numerical procedure works on blade row configurations of 
equal pitch. The actual SSME stator/rotor blade count Is 41:63, which is close 
to the 2:3 count shown In figure 8. A future goal of this work is to use this 
solution procedure to solve problems of unequal blade count like the one shown 
here. 


REFERENCES 

1. Chlma, R.V.: Development of an Explicit Multigrid Algorithm for Quasi 

Three-Dimensional Viscous Flows in Turbomachinery . AIAA Paper 86-0032, 
Jan. 1986. (Also, NASA TM- 971 28 . ) 

2. Jameson, A., et al.: Numerical Solutions of the Euler Equations by Finite 

Volume Methods Using Runge-Kutta lime-SteppIng Schemes. AIAA Paper 
81-1259, June 1981 . 

3. Sorenson, R.L.: A Computer Program to Generate Two-Dimensional Grids 

About Airfoils and Other Shapes by the Use of Poisson's Equation. NASA 
TM-81198, 1980. 
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Figure 7. - Static pressure 
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SIMULATION OF MULTISTAGE TURBINE FLOWS 

John J. Adamczyk 
NASA Lewis Research Center 
Cleveland, Ohio 

and 

Richard A. Mulac 
Sverdrup Technology, Inc. 
Mlddleburg Heights, Ohio 


A flow model has been developed for analyzing multistage turbomachinery 
flows. This model, referred to as the "average passage" flow model, describes 
the time-averaged flow field with a typical passage of a blade row embedded 
within a multistage configuration. The presentation summarizing the work done 
to date, based on this flow model, will be In two parts. The first part of the 
talk will address formulation, computer resource requirement, and supporting 
empirical modeling, and the second part will address code development with an 
emphasis on multitasking and storage. The presentation will conclude with 
Illustrations from simulations of the space shuttle main engine (SSME) fuel 
turbine performed to date. 
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POWERHEAD COMPONENT ARRANGEMENT 
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PREMISE 


OBSERVATION 


o HIGH SPEED MULTI - STAGE TURBOMACHINERY FLOWS 


HAVE TOO MANY LENGTH AND TIME SCALES TO 
BE AMENABLE TO DIRECT NUMERICAL SIMULATION 
EVEN ON TODAY'S MOST ADVANCED COMPUTERS 


o MODELS OF MULTI - STAGE FLOWS WHICH GIVE AN 
"AVERAGED" DESCRIPTION OF THE FLOW WITHIN 
TURBOMACHINERY PROVIDE USEFUL INFORMATION 

Figure 2. 


QUESTION 

GIVEN TODAY'S COMPUTER RESOURCES AND 
HIGH RESPONSE INSTRUMENTATION, IS IT TIME 
TO DEVELOPE MODELS WHICH PROVIDE A HIGHER 
DEGREE OF RESOLUTION OF MULTI - STAGE FLOWS 
THAN TODAY'S AXISYMMETRIC MODELS? 

Figure 4. 


MOST MODELS CURRENTLY USED TO ANALYZE 
MULTI - STAGE FLOWS ARE BASED ON AN AXI- 
SYMMETRIC REPRESENTATION OF THE FLOW 
WITHIN 1HESC MACHINES 

Figure 3. 


CONSTRAINT 


o PROPOSED MODEL MUST BE COMPATIBLE WITH 
THE AVAILABLE COMPUTER RESOURCES AND 
INSTRUMENTATION LIMITS 


o PROPOSED MODEL MUST HAVE A RATIONAL 
BASIS 


Figure 5. 
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TURBOMACHINERY MODELING EQUATIONS 


NAVIER 
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FLOW FIELD 


QUASI 
1-D EQNS. 


Figure 6. 

THREF— DIMENSIONAL CELL-CENTERED FINITE VOLUME 

FLOW SOLVER 


ADAMCZYK'S AVERAGE-PASSAGE EQUATION SYSTEM 
(d\u/dt) + L(\u) + f \£dVol - J XKdVol 
yj - (p ,pv r , rpv, , pv t , pe 0 ) 


L(Xu) - (AftL4 r + XG*dAj + XIMA^) 

*Lfl 

XSdVol - body forces, energy sources, momentum and energy 
J temporal correlations associated with neighboring 

blade row ( closure term) 


~ source term due to cylindrical coordinate system 
FOR ROTATING STSTEMS 

(dAw/<*0labs - (d\y/dt)\rel - Q(d\u/dO)\rel 
L(Xu) - I (XFhIA r + X(C - rQu)*dLA f <■ \H«1A # ) 


Figure 7. 
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CLOSURE STRATEGY 


FIELD EQUATION 


( a u 0) / dt)dLV + I ( u <0 ) + J S ( °d K = 0 
(du ( °/ 3i )dE + X(u a1 ) + 


5 W d K = 0 


Source Term 


Figure 9. 


CLOSURE STRATEGY 

BLADE ROW ( 1 ) 

X ( u 1 ^ ) + J S^dV - 0 

BLADE ROW ( 2 ) 

~L(u ( £) + J S (l) dV » 0 

WHERE 

S >,( l - 5 <0,t | Body Force, Energy Src, Velocity Cor., Energy Cor.) 


Figure 10. 
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ASSUME 
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UPON CONVERGENCE 
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Figure 11 . 


Figure 12. 


MODULAR CODE CONSTRUCTION 



Figure 13. 


) A dV - 0 
) A d V - 0 
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AVERAGE-PASSAGE EQUATION SYSTEM 



OUTPUT \ 


Figure 14. 


MULTITASKING OF MULTISTAGE 
3-D FLOW FIELD CALCULATION 




avtt xmp 



CPU II 
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Figure 15. 
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EVOLUTION OF THE TOTAL TEMPERATURE 
FIELD WITHIN THE S.S.M.E. FUEL TURBINE 



1st VANE 
MID PASSAGE 



1st ROTOR 
MID PASSAGE 



2nd VANE 
MID PASSAGE 



2nd ROTOR 
MID PASSAGE 


SIMULATION PERFORMED ON LEWIS CRAY XMP 24 


Figure 16 . 
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PROGRESS IN THE PREDICTION OF UNSTEADY HEAT TRANSFER 
ON TURBINES BLADES 
T. Cebecl 

California State University 
Long Beach, California 

R.J. Slmoneau 

NASA Lewis Research Center 
Cleveland, Ohio 

and 

A. Kralner and M.F. Platzer 
Naval Postgraduate School 
Monterey, California 


This presentation describes the progress toward developing a general 
method for predicting unsteady heat transfer on turbine blades subject to 
blade-passing frequencies and Reynolds numbers relevant to SSME, such as Illus- 
trated In figure 1. The method employs an Invlscld/vlscous Interactive proce- 
dure (fig. 2) which has been tested extensively for steady subsonic and 
transonic external airfoil problems. One such example Is shown In figure 3. 

The agreement with experimental data and with Navler-Stokes calculations yields 
confidence In the method. The present work extends the technique to account 
for wake generated unsteadiness. The flow reversals around the stagnation 
point caused by the nonuniform onset velocity are accounted for by using the 
Characteristic Box scheme developed by Cebecl and Stewartson. The coupling 
between the Invlscld and viscous methods Is achieved by using a special proce- 
dure, which, with a novel Inverse finite-difference boundary-layer method, 
allows the calculations to be performed for a wide range of flow conditions. 
Including separation. 

Preliminary results are presented for the stagnation region of turbine 
blades for both laminar and turbulent flows. A laminar model problem corre- 
sponding to a flow on a circular cylinder which experiences the periodic pass- 
ing of wakes from turbine blades Is presented to demonstrate the ability of the 
method to calculate flow reversals around the stagnation region. Calculations 
cover a Strouhal number range of 0.1 to 0.2 and a leading-edge Reynolds number 
range of 2xl0 4 to 4xl0 5 . 


BASIC EQUATIONS 

The basic equations are presented In figure 4. As an essential prelimi- 
nary to more extensive calculations Involving complete blades, we consider here 
the development of an unsteady boundary-layer method for calculating the flow 
properties near the stagnation region of a blade where the movement of the 
stagnation point with space and time poses problems. We assume that the exter- 
nal velocity distribution Is represented In dimensionless form by a function 
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(eq. (1)) which allows the variation of the stagnation point and the free 
stream velocity. 

For a two-dimensional, incompressible time- dependent laminar flow, the 
boundary layer equations and their boundary conditions are well known, and for 
conditions with no mass transfer and specified wall temperature they can be 
written in the form shown in equations (2) to (6). 

The calculation of upstream boundary conditions in the (t,y) plane at sortie 

x = x 0 when the conditions at a previous time line are known can introduce 

different problems. To illustrate these difficulties for the case of a moving 
stagnation point, let us consider equation (1). Since u e - 0 at the stag 
nation point by definition, its location based on the external streamlines is 
given by £ s = B(t) . 

Figure 5 shows the variation of the stagnation point with time according 
to the equation with B(t) = 1 + c sin ut, c = 1 , <j - ir/4 . We see that the 
stagnation point is at 2 when t = 2 and at 0 when t = fa, etc. 

If £ s were fixed, we could assume that u - 0 at Es = -1 for all time 

and for all y, but this is not the case. 

It is more convenient and useful to express equations (2) to (5) in a form 
more suitable for computation, lo achieve this, we introduce the dimensionless 
variables t, n, w, m, and G together with a dimensionless stream function 
f(S,n,T) and, with 0 - af/3^, write equations (2) to (5) as shown in equa 
tions (7) to (9). We use the Characteristic Box method to avoid the numerical 
problems associated with flow reversals in the stagnation region and to gener- 
ate the initial conditions on the next time "line." This scheme requires that 
equations (7) and (8) are expressed in terms of new coordinates. For this pur- 
pose, we note the definition of local streamlines and let dr = d^/f ' . If 
the distance in this direction is designated by s and the angle that it makes 
with the T-axis by 6, then the transformed momentum and energy equations (7) 
and (8) and their boundary conditions can be written as shown in equations (10) 
to (13). 


MODtl PROBLEM 

To Illustrate the method we introduce a model problem, which corresponds 
to a flow on a circular cylinder of diameter 0, which experiences the periodic 
passing of the wakes from the turbine blades (see fig. 6). For a time period 
tg the cylinder is subjected to a freestream velocity and for t w it is 

Immersed in a superimposed moving wake which has a motion component up. lhe 
cycle repeats Itself with a so-called blade passing frequency F[E T /( tg + t w ) ] 
and is related to the Strouhal number St by FD/U<». We assume that the stag- 
nation region of the cylinder is subjected to a velocity which varies in space 
and time according to equation (1). (The variables are listed in fig. 6.) 

Calculations are made for two values of Strouhal number taken equal to 0.1 
and 0.2 with the dimensionless minimum wake velocity u m = 1/3 in both 
cases. The computed values of wall heat flux G‘(0) show that they are not 
Influenced by the changes in the freestream velocity and are virtually constant 
for the range of E and t values considered with G 1 ( 0 ) ~ 0.50 for 
St = 0.1 and G'(0) ~ 0.51 for St = 0.2. 
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As shown In figure 7, the computed values of wall shear f“(0) for 
St = 0.1 are significantly Influenced by the changes In the freestream veloc- 
ity which causes flow reversals In the velocity profiles around the stagnation 
point based on the vanishing of the external velocity. The movement of the 
stagnation point and the resulting flow reversals Increase with time and with 
space. For example, the calculations for steady state have the stagnation 
point at 5=0, and, as expected, there Is no flow reversal In either side 
of the stagnation point. At t = 0.05, the stagnation point moves to 
5 = 0.15 but the flow reversals In the velocity profiles continue up to and 
Including 5 = 0.85 as can be seen from the results shown In figure 8. 


CONCLUDING REMARKS 

The computed results show that the numerical procedure Is able to obtain 
solutions for a range of blade-passing frequencies of practical relevance. The 
movement of the stagnation point with space and time and the resulting flow 
reversals around the stagnation point cause no computational difficulties and 
the numerical tests show that the accuracy Is better than required for practi- 
cal problems. 


GENERAL 0BJFCT 1VF 

TO DEVELOP A GENERAL METHOD FOR COMPUTING UNSTEADY HEAT TRANSFER ON TURBINES BLADES 
SUBJECT TO NONUNIFORM ONSET VELOCITY 



APPROACH 

INTERACTIVE BOUNDARY-LAYER THEORY BASED ON THE SOLUTION OF INVISCID AND VISCOUS FLOW 
EQUATIONS BY A NOVEL COUPLING PROCEDURE 


Figure 1. - Approach to rotor wake heat-transfer problem. 
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Figure 2. - Interactive scheme. Viscous 
calculations on the airfoil and in the 
wake are performed by a series of 
successive sweeps to establish new 
boundary conditions for the inviscid 
method . 
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Figure 3. Typical results for single 
airfoi Is. 
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GENERAL FORM OF EXTERNAL VELOCITY FIELO: 
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BOUNDARY-LAYER EQUATIONS FOR UNSTEADY FLOWS: 
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Figure 4. - Description of the method. 
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PRESENT CONTRIBUTION ADDRESSES THE NUMERICAL SOLUTION OF THE BOUNDARY-LAYER 
EQUATIONS FOR A GIVEN PRESSURE DISTRIBUTION AND ITS EVALUATION FOR MODEL LAMINAR AND 
TURBULENT FLOWS. 

FOR NONUNIFORM ONSET VELOCITIES. THE STAGNATION POINT LOCATION t, VARIES WITH SPACE 
AND TIME. I.E. FOR AN OSCILLATING AIRFOIL WHOSE u # NEAR THE LEADING EDGE IS 

2 + 2 

u ~ — (1 + c sin ox) 

6 JT77* 



Figure 5. - Variation of stagnation point with time. 



A(t) 
B(t ) 

h 


m 2 + LQ(st) 2 g - hi 2 1 1/2 

h 

tan' 1 [ 3 . 3St (i - 1)] 

1 + u 1 - u 

— 2 ♦ — 2 C06[tt(1 - 100 Stx)] 

Figure 6. - Rotor wake model. 
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EXPERIMENTAL AND THEORETICAL STUDIES OF I IME- AVERAGED AND 

1 1ME- RESOLVED ROTOR HEAT TRANSFER 

M.G. Dunn, W.J. Rae, and D.L. Rigby 
Calspan University of Buffalo Research Center 
Buffalo, New York 


CUBRC efforts In support of the SSME structural durability program have 
concentrated on obtaining detailed time averaged and time- resol ved (or phase- 
locked) measurements on a full-scale rotating turbine both with and without 
cold gas injection and on theoretical studies designed to improve the predic 
tion capability for these turbine flows. The experimental efforts have concen 
trated on use of the Garrett TFE 731-2 hp turbine. However, It has been 
possible to apply the theoretical efforts to predicting heat flux distributions 
obtained for two additional turbines - i.e., (1) the Garrett low aspect ratio 
turbine (LART) and (2) the Teledyne 702 turbine. In performing these calcula 
tions, the ground rules have been that only the turbine Inlet parameters and 
the turbine geometry can be altered In performing the calculations for the 
respective turbines. 


EXPERIMEN1 

The experimental technique is the short duration, shock tunnel approach, 
in which fast response, thln-film resistance thermometers are used to measure 
surface temperature histories at prescribed locations on the turbine component 
parts. Heat-flux values are then Inferred from the temperature histories by 
using standard data reduction procedures. Ihe turbine being used Is the 
Garrett TFE 731 2 hp stage, and both the nozzle guide vanes and the rotor 
blades are heavily Instrumented with thin film heat flux gauges. A detailed 
description of this Instrumentation was previously given (ref. 1) at the 198b 
Structural Durability meeting. Depending on how the data from a particular 
heat-flux gauge are recorded, one can get either time resolved (or phase 
locked) or time-averaged results. Both types of data are Illustrated in this 
presentation. 


PREDICT IONS 

It is Important to be able to predict both the time- resol ved and the time 
averaged heat- transfer distributions described previously. For the time 
averaged prediction, a wide variety of predictive techniques have been 
developed that range from simple flat-plate correlations to computer codes 
based on the full three dimensional Navier Stokes equations. Between these 
limits, there exists a substantial body of analytical and numerical methods 
which account for most of the dominant physical phenomena while still being 
sufficiently convenient for use by a large fraction of those Involved in the 
field. As part of this NASA effort, several of these midrange methods have 
been assembled to perform state-of the-art predictions for several different 
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turbines. The specific turbines for which these predictions have been per 
formed are (1) the Garrett TFF. 731-2, (2) the Garrett low aspect ratio turbine 
(LARI), and (3) the Teledyne 702- hp turbine. Results from the first two tur 
bines were reported in reference 2 and from the third In reference 3. lypical 
results are given in figures 1 and 2. 

Predicting the unsteady flow field is a significantly less developed area 
of research than is predicting the steady state flow field. Some effort Is 
ongoing within CUBRC to analyze and Interpret the unsteady thermal loads on a 
typical rotating turbine. At this time, the ability to perform phase locked 
heat flux measurements on a turbine rotor is more advanced than the ability to 
predict same. Figure 3 is a typical phase- locked heat flux distribution 
obtained at approximately 15 percent wetted distance on a turbine blade sue 
tlon surface for a rotation speed of 2/, 000 rpm. 
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INLET STANTON NUMBER 



Figure 1.- COMPARISON OF EXPERIMENTAL DATA WITH PREDICTION FOR ROTOR 
BLADE OF TFE 731-hp TURBINE . 



Figure 2. -STANTON NUMBER DISTRIBUTION FOR TELEDYNE 702 BLADE. 
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N87-22771 


ANALYTICAL AND EXPERIMENTAL STUDIES OF FLOW- INDUCED VIBRATION OF 

SSME COMPONENTS* 


S.S. Chen, J.A. Jendrzejczyk, and M.W. Wambsganss 
Argonne National Laboratory 
Argonne, Illinois 


Components of the SSMEs are subjected to a severe environment that 
Includes high-temperature, high-velocity flows. Such flows represent a source 
of energy that can Induce and sustain large-amplitude vibratory stresses and/or 
result In fluldelastlc Instabilities. Three components are already known to 
have experienced failures In evaluation tests as a result of flow-induced 
structural motion. These components Include the liquid-oxygen (LOX) posts, the 
fuel turbine bellows shield, and the Internal Inlet tee splitter vane. We have 
considered the dynamic behavior of each of these components with varying 
degrees of effort: (1) a theoretical and experimental study of LOX post vibra- 

tion excited by a fluid flow; (2) an assessment of the Internal Inlet tee 
splitter vane vibration (referred to as the "4000-Hz vibration problem"); and 
(3) a preliminary consideration of the bellows shield problem. This paper sum- 
marizes our effort In the resolution of flow-induced vibration problems associ- 
ated with the SSMEs. 

1. An Assessment of LOX Post Vibration: The LOX post vibration problem 

was assessed. Various excitation mechanisms Including vortex shedding, fluid- 
elastic Instability, LOX flow-induced Instability, acoustic resonance, and tur- 
bulent buffeting were considered. A scoping experiment was designed and 
completed; the experimental data agree well with analytical results (ref. 1). 

It was concluded that the potential excitation mechanisms for unshielded LOX 
posts are fluldelastlc Instability and turbulent buffeting. 

2. A Mathematical Model for LOX Posts In Crossflow: A general model, 

based on unsteady flow theory, for LOX post vibration was developed. The model 
Includes both motion-dependent fluid forces and fluid excitation forces 
(refs. 1 and 2). In the past, studies of the fluldelastlc Instability of a 
tube array In crossflow have been based on three different flow theories: 
quasi-static flow theory, quasi-steady flow theory, and unsteady flow theory. 
The quasi-static flow theory predicts fluid-stiffness -controlled Instability 
only and Is applicable for high reduced flow velocities. (Reduced flow veloc- 
ity Is equal to the flow velocity divided by oscillation frequency and cylinder 
diameter.) The quasi-steady flow theory predicts both fluid-stiffness- and 

f luld-damplng-controlled Instability. It Is applicable for high reduced flow 
velocities; at low flow velocities It does not predict some of the Instability 
characteristics. Only a complete unsteady flow theory Is applicable In all 
cases . 


*Work performed for NASA under an agreement with the U.S. Department of Energy. 
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3. Fluid Excitation Forces: The fluid excitation forces acting on a 

square array were measured as a function of Reynolds number, Incoming I low 
conditions, and tube location (ref. 3). The main conclusions are as iollows: 

(1) The fluid forces acting on the first few rows depend on the incoming 
flow conditions. Once the flow has passed through 3 to 4 rows, the I low and 
Its excitation on the tubes reach a steady state condition. 

(2) In a square array, the RMS values of the lift coefficients are larger 
than the drag coefficients. 

(3) As the flow passes through a tube array, both the turbulence intensity 
and the resulting fluid excitation forces Increase. 

(4) The power spectral densities of the fluid excitation forces are fairly 
flat for reduced frequencies of less than O.T; they reach a peak at the 
Strouhal frequency and then decrease drastically with an increase in reduced 
frequency. 

(5) The convection velocity for drag force is higher than that for lift 
force; It varies from 0.5 to 0.7 of the gap flow velocity. 

4. Motion Dependent Fluid Forces: Motion dependent fluid force components 

gj and hj acting on a group of cylinders In crossflow are given in figure 1. 
Components gj and hj consist of fluid Inertia, fluid damping, and fluid 
stiffness forces. F luid- damping and fluid stiffness coefficients for a tube 
row in crossflow were obtained as a function of reduced flow velocity (ref. 4). 
Typical results are shown in figure 2 as a function of the reduced flow veloc 

1 1 y . In addition, the effects of oscillation amplitude, tube alignment, and 
flow velocity were studied. From the characteristics of the motion dependent 
fluid forces, the stability criteria and the effect of various system param 
eters can be determined. 

5. Splitter Vane Vibration: A preliminary assessment of the excitation 

mechanism(s) associated with a 4000 Hz vibration of the Internal Inlet tee 
splitter vane was made from the available Information (unpublished data of 
Chen's, 1985). The most probable mechanism responsible for the 4000 Hz vibra 
tion is a "lock-in" oscillation Involving the splitter vanes and the vortex 
formation In the near-wake of the trailing edges of the vanes. A combination 
of structural and fluid-dynamic attenuation methods is recommended to eliminate 
the detrimental vibration. 

Even though there are significant gaps In our current knowledge, designers 
have been able to put together complex, Interactive components, such as those 
of the SSMEs, which operate In severe environments (high temperatures, pres 
sures, and flows) and still provide useful service without significant prob 
lems. In many cases, some of the obvious flow induced vibration effects can 
be avoided by using common sense and experience without detailed study and cal- 
culations. However, In the SSMEs, several vibration problems have been Identl 
fled which require further consideration. As an example, the shields on the 
liquid oxygen posts, which have served to reduce the flow Induced vibration 
problems of the posts to an acceptable level, have caused severe additional 
loads on other structural components. It would be desirable to redesign the 
Injectors with a detailed consideration of flow Induced vibration effects so 
that the flow shields could be removed. Clearly, there is a need to develop 
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and perform a prioritized research agenda In flow-induced vibration which will 
lead to a comprehensive technology base consisting of design data, prediction 
methods, and computer codes required In the design and evaluation of SSME 
components . 
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Figure 1. - Schematic of a group of circular cylinders In crossflow and fluid 
force and displacement components. 




Figure 2. - Fluid-damping coefficients for a row of cylinders obtained at 
several flow velocities denoted by different symbols. 
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OVERVIEW OF THE INSTRUMENTATION PROGRAM 


William C. Nleberdlng 
NASA Lewis Research Center 
Cleveland, Ohio 


As you have or will hear In the rest of this conference, one of the major 
goals of the durability program Is to develop computational models of the 
structural behavior, the fatigue life, and the environmental conditions perti- 
nent to advanced reusable space propulsion systems. In many cases, however, 
there Is very little experimental data available to either verify these models 
or to use as Input conditions to the models. 

The Instrumentation portion of this durability program Is thus aimed pri- 
marily at developing sensors and measurement systems capable of obtaining the 
data necessary for the verification of computational models. Indeed, some of 
these sensor systems may and probably will find their way Into the test bed 
engine and possibly even Into flight engines. The major effort, though, Is to 
develop what Is needed for bench mark experiments In research rigs to more 
fully define the actual conditions experienced In such an engine. In some 
cases, the research rig Itself will turn out to be the test bed engine simply 
because the conditions cannot be adequately simulated In less extensive 
facilities . 

One of the characteristics of measurement systems needed to verify codes 
Is that the sensors must be totally nonlntruslve or at least minimally Intru- 
sive so as not to significantly perturb the conditions being measured. The 
reason for this Is that the computer codes do not take the perturbations caused 
by something like a ppobe Into account. This leads to a heavy emphasis on 
laser optical techniques and on minimally Intrusive sensors. 

Another characteristic of such Instruments Is that they must be highly 
accurate and produce very high spatial and temporal resolution of the parameter 
being measured. Verifying codes demands this because the codes usually calcu- 
late highly resolved maps of such things as velocity, strain, temperature, etc. 

The measurement systems needed generally fall Into two broad categories. 
First there are the measurements needed on the surfaces of components such as 
turbine blades and vanes. Some of the desired parameters are temperature, 
strain, and heat flux. Papers on these three measurements are Included In this 
session. The other broad category encompasses those measurements needed In 
the flow environment around these components. Here, the desired results are 
highly resolved maps of such parameters as flow velocity, temperature, density, 
pressure, and species concentration. One of the papers In this session deals 
with this type of measurement. 

In summary, the Instrument program Is aimed at high-accuracy, high- 
resolution, minimally Intrusive measurements on space propulsion system compo- 
nents. The primary goal Is code verification, but some of the results may 
apply to later flight Instruments. 
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PROGRESS ON THIN-FILM SENSORS FOR SPACE PROPULSION TECHNOLOGY 

Walter S. Kim 
NASA Lewis Research Center 
Cleveland, Ohio 


The objective of this work Is to develop thln-fllm thermocouples for SSME 
components. Thln-fllm thermocouples have been developed for aircraft gas tur- 
bine engines and are In use for temperature measurement on turbine blades to 
1800 °F. The technology established for aircraft gas turbine engines will be 
adapted to the materials and environment encountered In the SSME. Specific 
goals are to expand the existing In-house thln-fllm sensor technology and to 
test the survivability and durability of thln-fllm sensors In the SSME 
environment. 


STATUS OF WORK 

An Initial goal of thln-fllm sensor work was to augment the In-house capa- 
bility In thln-fllm sensor technology. A thln-fllm sensor laboratory was 
established In a refurbished clean room which has been equipped with a new 
sputtering system for sputter deposition of both Insulator and sensor films 
(fig. 1). Other equipment Installed In the clean room Includes a photolitho- 
graphy exposure system for masking sensor patterns, parallel gap welder, and 
ultrasonic wedge bonder for attaching leadwlres to sputtered thin films. A 
scanning electron microscope was acquired for surface characterization of 
sputtered films and substrates. Also, a vacuum furnace was acquired to grow 
oxide coatings in controlled environments. 

One major accomplishment was the fabrication and testing of a thln-fllm 
thermocouple on SSME HPFTP blade material. The thln-fllm thermocouple was 
fabricated on a sample of the HPFTP blade material (MAR M-246 + Hf) following 
the fabrication process used for aircraft engine turbine blade materials. The 
test specimen was coated with an MCrAlY anticorrosion coating and oxidized to 
grow an adherent, electrically Insulating surface film of aluminum oxide. This 
film was augmented with additionally sputtered aluminum oxide upon which a 
type S thln-fllm thermocouple was sputter deposited (fig. 2). Junctions 
between the thin films and the leadwlres were formed by parallel gap welding. 

This thln-fllm thermocouple was tested for survivability In severe thermal 
shock conditions of 2000 °F and liquid nitrogen. This survivability test con- 
sisted of soaking the test piece at 2000 °F for 1 hour and then Immediately 
Immersing It In a pool of liquid nitrogen. The thln-fllm thermocouple survived 
five of these temperature excursions. After these Immersion tests, the sensor 
output signals were recorded while thermal shock cycling the sensor between 
2000 °F and liquid nitrogen (fig. 3). Sensor outputs were measured while con- 
secutively soaking and Immersing the sensor between those temperature extremes 
(fig. 4). The thln-fllm sensor has survived at least six of these thermal 
shock cycles. 
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Another major accomplishment was the fabrication of a thin film thermo 
couple on a future candidate SSMt material (fig. 5). I he thin- f 1 lm thermo 
couple was deposited on a flat specimen of a single- crystal nickel base 
superalloy such as PWA 1480. The preliminary tests Indicated that this sensor 
is just as good as the sensor prepared on MAR M 26 (t-Hf) specimen, lests are 
underway to determine the survivability and durability of this sensor in the 
same thermal shock cycling conditions. However, the severity of the thermal 
shocks Imposed on these sensors are not as severe as expected in SSMt and thus 
the next step is to test these sensors in the SSMt environment simulated in the 
Turbine Blade Tester at MSFC. 

The original plan for testing thin film sensors in an SSMt environment 
has been modified. The plan to contract for test hardware to test aircraft 
engine thin-film thermocouples was changed so that in house fabricated 
thermocouples on SSMt blades will be tested. The current plan is to install 
thin-film thermocouples on f 1 ight- hardware blades and other blades made from 
candidate material such as single-crystal nickel-base superalloy. We obtained 
from MSFC a blade holder assembly that was previously used for a thermal 
barrier coating experiment. Thin film sensors will be installed on this 
assembly and will be tested in the Turbine Blade Tester at MSFC. 


SUMMARY 

Significant progress was made in thin film sensor work by completing the 
establishment of a thin film sensor laboratory. One accomplishment was the 
development of a thin film thermocouple on SSMt HPFTF 1 blade material. Ihis 
sensor was thermal shock cycled for 12 cycles of consecutively soaking the test 
specimen between 2000 °F and liquid nitrogen. Another accomplishment was the 
development of a thin-film thermocouple on a candidate SSMt material. This 
thin-film thermocouple has been deposited on single crystal material of nickel 
base superalloy, and tests are underway to determine the survivability and 
durability under thermal shock cycle conditions. Lastly, the original plan to 
contract to test thin-film thermocouples on turbine blade material in an SSMt 
environment has been changed. The current plan is to install thin film thermo 
couples on f 1 Ight- hardware blades and single crystal blades mounted on the 
blade holder assembly obtained from MSFC for testing in their lurbine Blade 
Tester. 
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Figure 1. - Thin-film sensor laboratory. 
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Basic thin-film thermocouple technology. 
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HIGH TEMPERATURE AND PRESSURE LASER DIAGN0S1 ICS FOR SSME APPLICATIONS 


J.W.L. Lewis 

University of Tennessee Space Institute 
Tullahoma, Tennessee 


No paper available at time of printing. 
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PROGRESS IN OPTICAL STRAIN MEASUREMENT SYSTEM DEVELOPMI N I 


Christian T. Lant and Walld Qaqlsh 
Sverdrup Technology, Inc. 
Mlddleburg Heights, Ohio 


A laser speckle strain measurement system has been built and tested for 
the NASA Lewis Research Center. The system Is based on a speckle shift tech 
nlque, which automatl ca 1 1 y corrects for error due to rigid body motion, and 
provides a near real time measure of strain. This paper covers the first of a 
multiphase effort In developing an optical strain gauge capable of mapping In 
two dimensions the strain on the surface of a hot specimen. The objectives of 
this first phase have been to provide a noncontact, one dimensional, differen 
1 1 a 1 strain gauge for experimental purposes, and to determine the maximum open 
air temperature limit of the system. 


BASIC PRINCIPLES 

Laser speckle Is a phase effect caused by the diffuse reflection of spa- 
tially coherent light off of a rough surface. Its characteristically mottled 
Intensity distribution Is present In most laser applications. This optical 
strain measurement system Is based on the speckle shift technique of Ichlrou 
Yamaguchl (ref. 1), which utilizes the linear relationship between surface 
strain and the differential shift of speckle patterns generated by symmetri 
cally Incident laser beams. The speckle patterns created by the test specimen 
are Interpreted as higher order Interference fringes generated by a random dif- 
fraction grating, which Is the naturally rough surface of the specimen. Strain 
Induced on the specimen causes a change In spacing of the surface features, 
which In turn shifts the position of the Interference pattern (speckles). 

Theory shows that the displacement of a speckle pattern generated by a 
deformed specimen surface contains terms of translation, rotation, and strain, 
from which the strain term must be extracted. This extraction relies on the 
use of two symmetrically Incident laser beams reflected sequentially onto a 
linear photodiode array. After the reference (before strain) and shifted 
(after- strain) speckle patterns from each beam are correlated, taking the dif- 
ference In shift between them leaves only the component due to surface strain. 

The system has a gauge length of 5 millimeters given by the laser spot 
size, and a resolution of 16 microstrain determined by the sensor pitch of the 
photodiode array and the specimen- to-sensor separation. The theoretical strain 
measurement error Is 118 microstrain ±0.3 percent of the strain reading. 


RESULTS 

Uniaxial loads were applied to Inconel 600 test specimens by a testing 
machine. A number of data runs were performed with the system, confirming the 
expected response of the system as well as providing first-hand knowledge about 
sensitive parameters. Tests conducted during Phase I verified the stability 
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of speckle patterns generated by an oxidized specimen subject to strains over 
0.1 percent. Isothermal stress- strain plots as well as measurements of thermal 
strain showed good results up to temperatures of 450 °C. Stress- strain curves 
generated by the optical system were compared to data taken from conventional 
resistance strain gauges. Ihe modulus of elasticity for each set of data was 
calculated using a least squares fit and compared for each run. Agreement of 
the optical data was typically within 3 to 10 percent of the measured or hand 
book values of Young's modulus. 

The largest measurement error observed In the testing was due to Incom 
plete cancellation of the out-of-plane term of rigid body motion. The result- 
ing error In strain was caused by the radius of curvature of the laser beams. 

It was observed that the cause of the excessive out of plane motion encountered 
during testing was poorly machined specimen grips. However, by Increasing the 
radius of curvature of the laser beams, the error due to out-of-plane movement 
and other rigid body motions can be completely cancelled. 


CONCLUSIONS 

The results of the test runs show a good response of the system to surface 
strain. The residual error due to out-of-plane movement of the test specimen 
can be eliminated through modifications to the optical system. Hysteresis of 
the optical data over a load cycle Is shown to agree very well with resistance 
gauge data. Strain data Is obtainable at temperatures up to 450 °C In the 
present open-air configuration. Reduction of time-varying thermal gradients 
near the specimen, which cause decorrelation of the speckle patterns at higher 
temperatures, can Increase the temperature range of the system considerably. 


REf tRENCES 

1. Yamaguchl , I: A Laser speckle Strain Gauge. J. Phy. E., Scl. Instrum., 

vol. 14, no. 11, Nov. 1981, pp. 1270-1273. 
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HEAT FLUX CALIBRATION FACILITY CAPABLE OF SSME CONDI I IONS 

Curt H. Llebert 
NASA Lewis Research Center 
Cleveland, Ohio 


There Is a need to more thoroughly characterize the hostile space shuttle 
main engine (SSME) turbopump environment. It has been estimated that component 
surface heat flux In the hot-gas environment Is about 10 MW/m?, and this is 
about 50 times that encountered In aircraft engines. Also, material tempera 
ture transients can be as high as 1000 K In about 1 second. These transients 
can cause durability problems such as material cracking. Fleat flux sensors 
placed In the turbopump components can partially characterize this environment 
by measuring surface heat flux. These heat flux data can be used to verify 
analytical-stress, boundary layer, and heat transfer design models. 

Preliminary plans were discussed at the first SSME durability conference 
for designing and fabricating a new facility for the calibration and durability 
testing of prototype heat flux sensors for the SSME (ref. 1). This facility 
Is necessary for assessment of new heat flux gauge concepts needed in the hos- 
tile SSME turbopump environment. 


CALIBRATION FACILITY 

Figure T is a photograph of the calibration facility. A 100 kW Vortek 
arc-lamp and Anorad high-speed positioning system are shown positioned in the 
test cell. A Hewlett Packard 900-236 minicomputer Is located In a control room 
behind the control room window shown In figure 1. The Interrelation of cali 
bratlon facility components is schematically presented in figure 2. The mini 
computer Is used to control a microcomputer built Into the arc-lamp power 
module. The same minicomputer also controls a microcomputer located In the 
positioning system apparatus. Prototype heat flux sensors are placed at 
desired positions in the arc-lamp beam. Thermocouples attached to the beat 
flux sensor produce voltages which are scanned and read In a multiprogrammer 
and then stored In the computer. Calculation of heat flux is then performed 
using the stored voltages after they have been converted to temperature 
readings . 


ARC-LAMP ACCEPTANCE IES1S 

The Vortek arc- lamp system design is described In reference 1. Figure 3 
presents a comparison of heat flux measurements made in the Lewis facility and 
previously on an Identically designed arc lamp located in the Vortek facility. 
(Although the Lewis facility was not completed at this time, there were suffi 
cient apparatuses available to make the test.) The same water cooled heat flux 
sensor was used to perform both measurements. The agreement of heat flux data 
obtained with the two lamps (fig. 3) Is satisfactory and demonstrates an abil 
1 ty to repeat the heat flux within experimental error. Heat flux measurements 
and calibrations In the Lewis facility will also be compared with calibrations 
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of radiant heat flux sources at Case Western Reserve University and Pratt and 
Whitney. The comparisons are part of a grant with Case Western Reserve 
University to develop a credible basis for heat flux calibration at high heat 
flux loads. 


TRANSIENT HEAT FLUX MEASUREMENT AND ACCURACY 

It is not feasible to use commercial water-cooled heat flux gauges on SSME 
turbine surfaces because these gauges are too big and water cooling is not 
practical. Therefore, heat flux gauge design must be reinvestigated with par- 
ticular emphasis given to minimally intrusive instrumentation - possibly thin- 
film or miniature plug- type sensors. To aid in design, a transient heat flux 
prediction model for heat flow with variable thermal properties has been ana- 
lytically and experimentally developed. The model characteristics were studied 
using transient temperature data taken with heat flux gauges designed at Lewis 
Research enter for flux measurement in rocket engines in 1960 (ref. 2). This 
study is guiding the design of heat flux gauges for SSME and suggests that SSME 
heat fluxes can be measured with an accuracy of about ±1 to 20 percent. 


REFERENCES 

1. Liebert, C.H.: Heat Flux Sensor Calibrator. Structural Integrity and 

Durability of Reusable Space Propulsion Systems. NASA CP-2381, 1985, 
pp. 195-198. 

2. Liebert, Curt H.; Hatch, J.E.; and Grant, R.W.: Application of Various 

Techniques for Determining Local Heat Transfer Coefficients in a Rocket 
Engine From Transient Experimental Data. NASA TN D - 2 7 7 , 1960. 



Figure 1. - View of calibration facility. 
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NUMERICAL STUDY OF THE EFFECTS OF BOUNDARY CONDI I IONS ON THE MEASUREMENT 


AND CALIBRATION OF GARDON- T YPL HEAT FLUX SLNSORS 


M. Krane and A. Dybbs 
Case Western Reserve University 
Cleveland, Ohio 


To monitor the high-intensity heat flux conditions that occur in the space 
shuttle main engine (SSME), It Is necessary to use specifically designed heat 
flux sensors. These sensors, which are of the Gardon type, are exposed on the 
measuring face to high- intensity radiative and convective heat fluxes and on 
the other face to convective cooling. To Improve the calibration and measure- 
ment accuracy of these gauges, we are studying the effect that the thermal 
boundary conditions have on gauge performance. In particular, we are studying 
how convective cooling effects the field inside the sensor and the measured 
heat flux. 

The first phase of this study involves a numerical study of these effects. 
Subsequent phases will involve experimental verification. 

A computer model of the heat transfer around a Gardon type heat flux sen- 
sor was developed. The results of this study will allow us to make corrections 
both In calibration and applications so that the input (hot-side) heat flux can 
be determined more accurately. 

Two specific geometries are being considered as shown In figure 1: 

(1) Heat flux sensor mounted on a flat-plate 

(2) Heat flux sensor mounted at the stagnation point of a circular 
cylinder 

Both of these configurations are representative of the use of heat flux sensors 
In the components of the SSME. 

The purpose of the analysis Is to obtain a temperature distribution as a 
function of the boundary conditions. Previous analyses of Gardon sensors have 
concentrated on the thin- foil part of the sensor, considering It to be attached 
to an Infinite heat sink. This Is a reasonable assumption at low temperatures, 
but at the high temperatures found In the SSME, this assumption Is question- 
able. The results of this analysis will determine (1) how the cooling of a 
sensor In such extreme environments approximates the "thln-foll" assumption and 
(2) how nonuniformities in the cooling (caused by convection cooling coeffl 
clents being space dependent) effect the measurement of the hot side heat flux. 

The equations for the general case are shown In figure 2. We have not 
Included in this formulation the thermal conductivity variations caused by the 
thermocouple wires or channels. 
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As a first test case, the following simplifications were made: 

(1) Constant radiant hot-side heat flux 

(2) Constant convection cooling coefficients and ambient temperatures 

(3) Symmetrical k- dl s trlbutlon (l.e., no thermocouple wire or channels) 

(4) Constant thermal conductivity 

The first assumption eliminates natural convection or re-radiation cool 
ing. We know how to handle both of these effects. Ihe second and third 
assumptions eliminate any o dependence In the problem so that It Is axlsym 
metric. Finally, the fourth assumption makes the problem a linear one. 

A finite-difference model was used to solve the equations shown in 
figure 2. Several cases have been studied. Ihe results are shown In figures 3 
to 5. Figure 3 shows the calibration curves for a typical sensor. Ihe top of 
the sensor has a convective heat transf er coefficient h (- 0 p - 10 KW/m2°K. 

This figure shows the results of varying the cooling or bottom side heat 
transfer coefficient. For these cases, the sensitivity of the sensor decreases 
with an Increase in the cooling heat transfer coefficient. This result is also 
shown in figure 4 where the sensor sensitivity is plotted versus the cooling 
heat transfer coefficient. 

The cases studied in figures 3 and 4 were for radiation and convective 
heat fluxes at the top face of the sensors to be of the same order of magni 
tude. If the convective heat flux Is small In comparison to the radiation heat 
flux, then the cooling heat transfer has little effect on the calibration 
curve. This is shown in figure 5 where the cases of figures 3 and 4 are com 
pared to the case of ht op = 0 KW//m 2 °K with varying cooling h's. 
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Two Geometries: 


1. FLAT PLATE 



F igure 1 . 


Heat flux sensor configurations to be considered 



GENERAL CASE 



B.C . : TOP : kj |I (r,6,0,t) - aQ i - h CT (T(r,e,0,t) - T CT ) 


SIDE: 


BOTTOM: 


CENTER: 


9T 


-a a (T 4 (r,0, O.t) - T 4 T ) 
(R .6,z,t) ="h rD (T(R ,0,z ,t) - T-p) 


k l 3r "CR' , ' , V U ’^ CR 

[ k 1 |I (r,6,H,t) =-h CB (T(r,0,H,t) - T C g) ; Rf rsR fl 

k, (r,0, H,t) =-h rB (T ,r,0,H ,t) - T rR ; 05 r« R, 


2 3z v ' ,v ’ "CB 

T(O,0,z,t) | < « 


I.C. 


T(r,0,z,O) = T CI (r,0,z) 


Figure 2. - General case and formulation. 
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Figure 3 . - Calibration curves for several values of h^. Note that 
hton = 10 KW/m 2 °K Is constant and uniform for all cases. 
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Figure 4. - Variation of sensor sensitivity with bottom- Note that 
^top = 10 KW/m2°K was fixed and uniform In all cases. 
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x top rad CONV j^2j 

Figure 5. - Comparison of calibration curves for (1) large top-side cooling 
solid curves (h^op * 10 KW/m 2 °K) - 5 curves and (2) small top-side 
cooling broken line (h^ 0 p = 0 KW/m2°K) - 5 separate curves all on same line. 


OVERVIEW OF FATIGUE, FRAC1URE, ANO CONSlllUllVE MODELING PROGRAMS Al LEWIS 


Gary R. Halford 
NASA Lewis Research Center 
Cleveland, Ohio 


The objective of the NASA Lewis Research Center's Structural Integrity and 
Durability program is to develop and verify cyclic constitutive stress strain 
and crack initiation life prediction models for materials of use in reusable 
space propulsion systems. Emphasis is on application to those hot gas path 
components that have demonstrated limited durability (i.e., turbopump blades, 
LOX posts, and the main combustion chamber liner). 

Our approach to meeting the rigorous demands of space propulsion systems 
has been to borrow from technology developed under aeronautical gas turbine 
sponsorship. Where the aerotechnology can be used directly, it is being 
applied; where modifications and extensions are required, these are being pur- 
sued; and where no technology exists, it is being developed. A revealing con 
trast between the severity of aeronautical gas turbine engine operation and 
that in the turbopumps of the space shuttle main engine (SSME) is Illustrated 
in figure 1 . 

The program is in its fourth year, and significant progress has been made. 
Figure 2 lists several major accomplishments that range from development of an 
understanding of the complex micromechanisms of cyclic deformation in advanced, 
directional, nickel base superalloys at high temperatures, to a new generalized 
mean stress model for fatigue, and to an accurate, shortcut structural analysis 
computer code for thermally loaded structures. 

Durability research programs currently in progress or planned for the 
immediate future are listed in figure 3. Results of some of these are being 
presented in the papers that follow in this session. Finally, programs planned 
tentatively for the more distant future are presented in figure 4 . 

Considerable interest has been devoted to the single crystal nickel base 
superalloy, PWA 1480, in the anticipation that it will provide superior proper- 
ties to the bi 1 1 -of-mater i a 1 , MAR M 246 c Hf, for use in turbopump blading. A 
consensus of much of the work to date on this alloy points up the substantial 
benefit of reducing voids. Reductions both in numbers of voids and their size 
should significantly Improve high-cycle fatigue resistance and aid in the alle- 
viation of hydrogen interaction problems. Progress toward this end is reported 
in this session by researchers from Rocketdyne and Carnegie-Mel Ion University. 
Their two programs, while not funded currently out of the Lewis Durability pro 
gram, certainly have a direct bearing upon our research efforts. 
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THERMAL GRADIENTS CYCLIC STRESS-STRAIN RESPONSE 



CRACK GROWTH REGIME8 



A*,. STRAIN INTENSITY 


Figure 1. - Thermal loadings, structural response, and durability. (Space 
shuttle main engine (SSME) compared with aeronautical gas turbine engine) 


Hydrogen Interaction effects on mechanical strength 

• Fatigue loading magnifies extent of hydrogen 
Interaction compared with tensile loading. 

• Important micromechanism of hydrogen Interaction 
postulated, l.e.. Internal voids act as traps 
for hydrogen and enhance embrittlement. 

High mean stress In high cycle fatigue (HCF) 

• Factor of two potential Increase In HCF strength 
by reducing void size by factor of four. 

Micromechanisms of cyclic deformation In turbopump blading alloys 

• Deformation mechanisms leading to low ductility at 
temperatures where H 2 Interaction Is greatest have 
been documented and understood. Improved ductility 
could lead to reduced H 2 Interaction and to Improved 
high cycle fatigue (HCF) resistance. 

Cumulative damage rules for HCF/low cycle fatigue (LCF) Interactions 

• HCF/LCF results for 316 SS, Haynes 188, PWA 1480 show 
substantial life-degrading Interactions 

• Simple cumulative damage Interaction equations verified 
and applied In HCF/LCF life extension tradeoff studies 

lhermomechanlcal fatigue (TMF) and relation to Isothermal fatigue 

• Isothermal fatigue has known limitations for 
predicting thermomechanical fatigue (IMF) lives 

• Newly proposed blthermal fatigue testing technique 
aids In simplifying IMF Investigations 

• Newly proposed TMF life prediction model based on 
application of total strain version of strain range 
partitioning and bithermal cycling 

Advanced structural analysis demonstrations 

• Successfully applied advanced constitutive models 
to analyses of complex mission cycles of SSME 
turbopump blading, both MAR-M 246 <■ Hf and PWA 1480 

• Verified accuracy of shortcut structural analysis 
procedure, 10 000 to 1 savings In CPU time 


Figure 2. Accomplishments of structural durability program. 
( Fatlgue/fracture/constl tutl ve modeling) 
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HCF/LCF Interaction under mission-related loading 

• Interspersed HCF and LCF loadings 

• TMF Interaction with HCF loadings 

• Generalized creep-fatigue Interactions 

Application of mechanistically-based constitutive and life models 

• Integration of aero-and SSME-based program results 

• Pick up where HOST has left off 

Macroscopic effects of high pressure H 2 on constitutive and life models 

• Competitive procurement package In preparation 

• Establish validity of structural durability models 
developed under program 

Structural analyses demonstration program 

• GE contract to apply most sophisticated models 
available to analyze three SSME components, 

Including blade and LOX post 

Thermal ratcheting damage analyses 

• Interaction of ratchet strains with TMF and HCF 

• Application to main combustion chamber materials 

Figure 3. - Current and near-term programs. 

( Fatlgue/fracture/constltutlve modeling) 
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Cyclic crack growth under thermomechanical loading conditions 

• Application of path Independent integrals to 
describe crack driving function; pick up where 
HOST has left off 

Application of structural durability analyses to test bed components 

• Improved input information from test-bed instrumentation 

to provide for more accurate structural durability assessment 

Cyclic constitutive and life model development for advanced blade concepts 

• Metal matrix composite (MMC) blade materials evaluation 


Figure 4. - Long term programs. 
(Fatigue/f rac tur e/cons ti tutive model mg) 
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DEF0RMA1 ION AND FATIGUE BEHAVIOR OF SSME 1URB0PUMP BLADE MATERIALS 


Walter W. Milligan and Stephen D. Antolovlch 
Georgia Institute of Technology 
Atlanta, Georgia 


Directionally solidified and single crystal superalloys which are intended 
for use as turbopump blade materials are anisotropic, both elastically and 
plastically. Therefore, isotropic constitutive models must be modified. Sev 
eral models which are now being developed are based on meta 1 lurglcal theories 
of deformation In these types of alloys. However, these theories have not been 
fully justified, and the temperature and strain regimes over which they may be 
valid are poorly defined. The objective of this work Is to study the deforma 
tlon behavior of the alloys. In order to determine the validity of these models 
and to thereby support the ongoing research efforts in solid mechanics. 


REFERENCES 

1. Milligan, W.W.; and Antolovlch, S.D.: Yielding and Deformation Behavior of 

the Single Crystal Superalloy PWA 1080. Metall. Trans. A., vol. 18, no. 1, 
Jan. 1987, pp. 85-95. 

2. Milligan, W.W., Jr.: Yielding and Deformation Behavior of the Single 

Crystal Nickel-Base Superalloy PWA 1480. (USAAVSCOM I R86-C- 18, Georgia 
Institute of Technology; NASA Grant NAG3 503) NASA CR 175100, 1986. 

3. Huron, E.S.: High Temperature Monotonic and Cyclic Deformation in a 

Directionally Solidified Nickel-Base Superalloy. (USAAVSCOM 1R 86 C 19, 
Georgia Institute of Technology; NASA Grant NAG3 503) NASA CR- 175101, 1986. 
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Figure 1. - PWA 1480 initial microstructure. 
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Classical Plasticity Theory 





Figure 2. - Constitutive modeling. 



Theory, Experiment Coincide Theory 

Experiment 


figure 3. - PWA 1480 at 1600 °F - after K.P. Walker and t. Jordan. 
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Figure 4. - Unified crystallographic approach. 

Crystallographic approach requires physical validation 

-Active slip systems 

-Cube slip only observed in y' 

-Slip planarity 
-Deformation mechanisms 

Figure 5. 
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must be "overcome" for substantial plasticity 


Figure 6. Deformation mechanisms. 


Possibilities: 

1. Looping 

2. Shearing 

3. Climb-assisted Bypass 
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Figure 8. - PWA 1480 yield strength (0.05 percent offset) versus 
temperature. 



Figure 9. - PWA 1480 modulus normalized 
yield; strength versus Inverse tem- 
perature. 
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Figure 11. - Shear offsets parallel to b. 
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Figure 12. - Typical tensile fracture surface. 
(D.S. Mar- M 246 tested at room temperature.) 
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Figure 13. - Low temperature 
(Copley- Kear, Huther- Repplch) . 
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Figure 14. - Y 1 Bypass at high temperatures 
(0.3 percent e p ) . 


Figure 15. - Interfacial network of disloca- 
tions developed at 1093 °C after about 
0.25 percent e p . 
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Figure 16. - High dislocation density at failure 
at 871 °C. 


Climb is rate-limiting step 

-Controlled by diffusion 


Figure 17. - High temperatures (Brown 
and Ham, bypass model). 
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0.5 percent/min 



50 percent/min 


Figure 20. - 927 °C, Ae p = 0.1 percent, N = 4. 
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Figure 21. - Mar-M 246 deformation 
structures after low cycle fatigue at 
927 °C. 


Original unavailable at time of printing. 
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Figure 22 . - PWA 1480; 20 °C, 
Aep = 0.05 percent, N = 4. 
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Figure 23. - Surface slip traces 
after low cycle fatigue at 
705 °C, Mar-M 246. 


Original unavailable at time of printing, 
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Deformation during creep, 

uniaxial tension and LCF essentially similar 

-By-pass at 927 °C 
-Shearing at Lower Temps 

Figure 24. 


Implications 

-Unified models justified 

-No shearing in (001) crystals at 927 °C 

-(123) tests to come 


Figure 25. 
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FATIGUE DAMAGE INTERACT ION BEHAVIOR OF PWA 1480 


Michael A. McGaw 
NASA Lewis Research Center 
Cleveland, Ohio 


The fatigue damage interaction behavior of PWA 1480 single crystal alloy 
has been experimentally established for the two-level loading case in which a 
block of low-cycle fatigue loading is followed by high cycle fatigue loading 
to failure. A relative life ratio N1/N2 (where N1 and N2 are the low 
and high-cycle fatigue baseline lives, respectively) of approximately 0.002 was 
explored to assess the interaction behavior. The experimental results thus far 
show evidence of a loading order interaction effect to a similar degree of 
detriment as has been observed in polycrystalline materials. Current genera- 
tion single crystal alloys in general, and PWA 1480 In particular, contain 
pores; Indeed, it was observed In all cases that specimen failure initiated 
from pores connected with or Immediately below the surface. Detailed fracto- 
graphlc and metallographlc studies are currently being made to assess the 
nature of the porosity in terms of Its effect on fatigue life. 


EXPERIMENTAL RESULTS 

Specimens of PWA 1480 were fabricated from slab material. The specimen 
geometry was solid cylindrical, having a l-1nch uniform axis gauge section and 
threaded ends. The specimens were within 10° of the 001 crystallographic 
direction . 

Low-cycle fatigue tests were conducted In strain control using a special- 
ized electrohydraul 1c tension compression testing system (fig. 1, ref. 1). The 
specimens were tested at room temperature using a computer generated sinusoidal 
waveform command of constant frequency (0.2 Hz). Control and data acquisition 
were performed by computer using a software system described In reference 2. 

The strain gauge used for strain measurement and feedback had a gauge length of 
0.24 Inch. High-cycle fatigue tests were conducted with the same equipment, 
although under load control. The test waveform was sinusoidal; the frequency 
was 300 Hz. 

The results of this baseline fatigue characterization testing are dis- 
played In figures 2 and 3. Failure In all tests Is defined as specimen frac- 
ture Into two halves. 

The Interaction tests were based on reference low-cycle and high- cycle 
fatigue life levels of 3000 and 1.5 million cycles-to-fallure, yielding a 
relative life level ratio of approximately 0.002. As the purpose in this pro- 
gram is to Identify potentially deleterious life behavior on a cumulative 
fatigue basis, the tests were conducted so that the low cycle fatigue loading 
was applied first, and the remainder of the test was carried out under high- 
cycle fatigue loading conditions until failure occurred, whereupon the remain 
Ing life fraction was measured. 
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The results of the Interaction testing are displayed In figure 1 (where 
the data are Interpreted In terms of a stress based life criterion) and 
figure 5 (where the data are Interpreted In terms of a total strain based life 
criterion). The solid curves In both figures are predictions made using the 
Damage Curve Approach (ref. 3), where the average of the low-cycle fatigue 
lives predicted for the applied cycle parameter (stress or total strain) was 
used for the low-cycle fatigue reference life N-| and the same procedure was 
used for obtaining the high cycle fatigue reference life N 2 . 


DISCUSSION AND CONCLUSIONS 

A large amount of scatter In the baseline fatigue behavior Is observed in 
figure 3, especially In the 10 6 to 10 8 life range. A number of experiments 
were performed at nominally the same completely reversed stress range; three 
failures were obtained In the 10 6 range while two "runouts" were observed at 
the middle 10 7 range. Scatter, although apparently not to the same degree, 
was observed In the low cycle fatigue range as well. A preliminary fracto 
graphlc analysis of the failed specimens revealed cracks Initiating at surface 
or near-surface porosity In all cases. The nature of the porosity varied 
widely In terms of pore size, geometry, and orientation with respect to the 
loading (and crystal growth) direction. Preliminary analysis also indicated a 
change In fracture surface appearance as a function of life: while the failure 

always Initiated at a pore In the case of low-cycle fatigue, the early propaga 
tlon was on two intersecting crystallographic planes; In the case of high-cycle 
fatigue the early propagation was on a singularly dominant crystallographic 
plane. The actual crack propagation phase of failure observed In these experi- 
ments Is for all practical matters zero In the context of cyclic life. The 
stress strain response observed at the onset of failure gave little or no Indi- 
cation of the Impending fracture. 

The Interaction experiments show an Interaction effect produced by cumu- 
lative loading, although the precise nature of the Interaction effect Is not 
well correlated by current cumulative damage rules, such as the Damage Curve 
Approach. A difficulty In the study of cumulative fatigue damage has been the 
compounded scatter observed In such studies; that Is, If the reference life 
levels are known to be within a factor of 2, the loading Interaction cannot be 
known to be anything less than a factor of 2. In fact, the factors of 2 com 
pound one another. Preliminary f ractographlc analyses reveal fracture surfaces 
which appear quite planar, being perpendicular to the loading axis for applied 
low-cycle fatigue ratios less than approximately 0.5. 

The f ractographlc appearance Is reminiscent of polycrystalline materials 
falling In high-cycle fatigue. This fracture morphology was never observed In 
the baseline experiments. For applied low-cycle fatigue life ratios of approx- 
imately 0.5 and greater, the fracture takes on a character similar to that 
observed In the baseline low cycle fatigue experiments. 

Detailed f ractographlc and metal lographlc studies are being made on the 
failed specimens In an effort to explain the Influence of pore size, geometry, 
and orientation on fatigue life, as well as Interaction behavior. 
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SATEUJTE COMPUTER SYSTEM 


Figure 1. - High-cycle fatigue/low-cycle fatigue materials test system. 
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Figure 4. - Results of Interaction testing. Data 
Interpreted In terms of stress-based life 
criterion. 
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Figure 5. - Results of Interaction testing. Data 
Interpreted In terms of total strain-based life 
criterion. 
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A NEW FORMULATION OF MEAN STRESS EFFECTS IN FATIGUE 

S.S. Manson 

Case Western Reserve University 
Cleveland, Ohio 

and 


K.R. Heldmann 
Steelcase Corporation 
Grand Rapids, Michigan 


A common method of treating the mean stress effect on fatigue life Is to 
displace the elastic line on a Manson-Coff In-Basquln diagram while retaining 
the position of the plastic line. Manson and Halford pointed out that this 
procedure Implies that mean stress significantly affects the cyclic stress- 
strain curve (ref. 1). Actually, however, they showed experimentally and by 
more general reasoning, that mean stress has little. If any, effect on the 
cyclic stress-strain curve. Thus, they concluded that It Is necessary to dis- 
place the plastic line as well as the elastic line In order to keep the cyclic 
stress-strain curve unaltered. Another way of expressing the common displace- 
ment of the two lines Is to keep the lines In place and change the horizontal 
coordinate to Include a term relating to the displacement. Thus, Instead of 
life, 2Nf, as the horizontal coordinate, a new coordinate can become 2Nf 
[1 - o m /of] 1/b , thereby displacing both the elastic and plastic lines by an 
amount [1 - o m /of]Ub where a m Is the mean stress and a'f Is the Intercept 
of the elastic line at N f = 1/2 cycles and b Is the slope of the elastic 
line. 


Such a procedure does not, however, always produce the proper effect on 
the Goodman diagram In which mean stress Is plotted against alternating stress 
for selected values of life. The Goodman diagram for this case would necessar- 
ily be straight lines for all materials. While some materials do, Indeed, pro- 
duce straight lines, others are known to be concave, while still others have 
convex curvature. In fact, for some materials, the Goodman diagram may be con- 
cave for one life level, straight for another life level, and convex for a 
third life level. A more generalized behavior can, therefore, be achieved by 
using, as the horizontal coordinate, the parameter 

/<j\A + B log N f l 1/b 
1 / _JB\ t 

Depending on the values of A and B, the curvature of the Goodman lines can 
then change according to life level. 
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We have examined a number of materials for which various degrees of curva- 
ture have been observed, as well as one material for which the curvature 
changed from convex to concave as life was changed. 

Figure 1 shows schematically the concept behind the procedure we have 
developed. Figure 1(a) shows the coordinate axes In a Manson-Cof f In-Basquln 
plot while figure 1(b) shows an example whereby curvature Is altered at dif- 
ferent values of life level for particular parameter values of A and B, as 

shown. 

Figure 2 shows the application of this type of analysis to various alumi- 
num and steel alloys for which extensive data are available In the literature. 
The lines drawn In these figures are derived from the diagrams of figure '(a) 
for the specific constants shown for each material. The experimental points 
are derived from faired curves through the data. Good agreement Is seen In all 
cases. Note especially the 9N1-4CO-0.45C material for which different curva- 
tures appear at different life levels. 


REFERENCE 


1. Halford, G.R.; and Manson, S.S.: Discussion of Paper: M. Doner, 

K.R. Bain, and J.H. Adams; Evaluation of Methods for the Treatment 
Stress Effects on Low-Cycle Fatigue. J. Eng. Power, vol . 104, no. 
1982, pp. 403-411: J. Eng. Power, vol. 104, no. 2, Apr. 1982, p. 


of Mean 
2, Apr. 
411 . 


90 



Figure 1. - Framework for new method 
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Figure 2. - Comparisons between observed and calculated axial fatigue strengths. 
( Predictions ; From experimental curves) 
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Figure 2. - Continued 
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NONLINEAR HEAT TRANSFER AND STRUCTURAL ANALYSES OF SSME TURBINE BLADES 

A. Abdul-Azlz 
Sverdrup Technology, Inc. 

Mlddleberg Heights, Ohio 

and 

A. Kaufman (Retired) 

NASA Lewis Research Center 
Cleveland, Ohio 


Three-dimensional nonlinear finite-element heat transfer and structural 
analyses were performed for the first stage high-pressure fuel turbopump blade 
of the space shuttle main engine (SSME). Directionally solidified (DS) 

MAR-M 246 material properties were considered for the analyses. Analytical 
conditions were based on a typical test stand engine cycle. Blade temperature 
and stress-strain histories were calculated using MARC finite-element computer 
code (ref. 1). This study was undertaken to assess the structural response of 
an SSME turbine blade and to gain greater understanding of blade damage mecha- 
nisms, convective cooling effects, and the thermal-mechanical effects. 

Hot-gas path components for reusable space propulsion systems operate 
under extreme gas pressure and temperature. These operating conditions subject 
the high pressure stage turbine nozzles and blades to severe thermal transients 
that can result In large Inelastic strains and rapid crack Initiation. 

Advances In casting techniques have allowed the development of directionally 
solidified and single crystal alloys for high temperature components In space 
propulsion vehicles. Mechanical anisotropy exhibited by these alloys have to 
be taken Into account In the analytical studies. To Improve the durability and 
accuracy of turbine blades and other hot section components, an accurate knowl- 
edge of the temperature and stress-strain histories at the critical location 
for crack Initiation Is required (fig. 1). 

Experimental measurements of gas temperature profiles In the SSME turbo- 
pumps are difficult to obtain because of the high gas temperature and severe 
thermal transients. Turbine blade temperatures are primarily a function of hot 
gas flow and cooling. The temperature field Is determined by the heat transfer 
from the hot gas to the blade. This heat transfer and Its variations are 
determined by knowledge of the gas film coefficients. Also, the time- 
temperature history profile at the start transient, steady state, and cutoff 
Is obtained through a combination of analytical and experimental results which 
Is due to the complex flow phenomena through an accelerating turbine. 

Temperature-dependent properties for the MAR-M 246 + Hf alloy were mainly 
provided by Rockwell International Corporation. These elastic properties are 
summarized In table I. Mean thermal coefficient of expansion data were con 
verted to Instantaneous values for MARC Input. Longitudinal stress-strain 
properties, summarized In table II, were used for elastic-plastic region. The 
mission used for this analysis Is shown In figures 2 and 3 In terms of Inlet 
temperature, gas pressure, and revolutions per minute (RPM). This cycle Is 
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applicable to a factory test of the engine; It Is also reasonably representa- 
tive of a flight mission except for the foreshortened steady-state operating 
time. The major factor Inducing fatigue cracking Is the transient thermal 
stresses caused by the sharp Ignition and shutoff transients. The finite- 
element model of the blade In figure 4 was constructed of a three-dimensional 
eight-node Isoparametric brick element. The model consisted of 1-25 elements 
with 1575 nodes and 4660 unsuppressed degrees of freedom. 

Heat transfer coefficients at the blade airfoil (figs. 5(a) and (b)) were 
predicted by running a boundary layer analysis using a modified version of the 
STAN5 boundary layer code (refs. 2 and 3). Thermal environment experienced by 
the platform and shank were obtained from reference 4. Details regarding pre- 
diction of heat transfer coefficients at the stagnation region for the airfoil 
are available In reference 5. Transient results were obtained by scaling 
steady-state heat transfer coefficients based on transient flow and tempera- 
ture. Predicted high temperature locations were evaluated for the two tempera- 
ture spikes shown In figure 2. 

The thermal response predicted from the finite-element analysis showed 
that the leading and trailing edges of the airfoil base are the hottest loca- 
tions. Temperature distributions showed a cordwlse variation at the first 
Ignition spike and a spanwlse variation thereafter Into the cycle. A uniform 
temperature distribution was dominant at most of the airfoil surface during 
cruise except near the base at the platform junction where a mixture of cold 
and hot gas Is present. The coldest spot was always at the blade root because 
of cooler boundary conditions. Blade temperature Is shown In figures 6 
through 8. 

Elastic-plastic analyses have been conducted for the HPFTB blade with MARC 
code. Plastic strain calculations were based on Incremental plasticity theory 
using Von Mlses Yield criterion, the normality rule and a kinematic hardening 
model. The material elastic-plastic behavior was specified by the yield 
strengths and work hardening properties In the longitudinal direction; trans- 
verse properties were not available. Creep analyses were not performed at the 
present time because of Inadequate knowledge of the creep characteristics for 
anisotropic blade material. 

Incremental loading Included centrifugal and gas pressure loads and metal 
temperature distributions as calculated from the heat transfer analysis. The 
same Increments were used for the heat transfer and the elastic-plastic struc- 
tural analysis. Approximately two million words of core storage on the CRAY- 
XMP computer were needed to run the problem. Analysis required about 4 hours 
of central processor unit (CPU) time on the CRAY system. The directionality 
of the elastic material properties causes anisotropic constraints. Lekhnltskll 
(ref. 6) has derived the generalized elastic strain equations for an aniso- 
tropic body with a transverse plane of Isotropy. This anisotropic stress- 
strain law was Incorporated In MARC user subroutine HOOKLW. 

The colder airfoil base temperatures Induce tensile thermal stresses at 
the critical leading edge location that are additive to the centrifugal 
stresses. A noticeable discrepancy Is seen In the compressive strain region 
where the reduction In elastic analysis points resulted In failure to capture 
some of the cycle fluctuation due to transient thermal effects during the rapid 
engine cooldown. The CPU time for the elastic finite-element analyses amounted 
to 5 percent of that required for one cycle of the nonlinear finite-element 
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analysis. The critical location for crack Initiation Is at the leading edge 
near the base. Results of both elastic and Inelastic structural analyses Indi- 
cated the region of the finite-element model with the largest total strain 
range was coincident with the observed crack Initiation site. 

Radial stress distributions as calculated at the gausslan Integration 
points closest to the suction and pressure surfaces are presented In figures 9 
through 11. The stresses primarily reflect the centrifugal and thermal load- 
ings. It Is apparent that the analysis would not show a low-cycle fatigue 
problem If transient thermal effects were not considered. 
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TABLE I. - OS MAR-M 2*6 PHYSICAL PROPERTIES 


Temperature, 

•c 

Modulus of elasticity, 
GPa 

Means coefficient of 
thermal expansion, 
percent/*C 

Longitudinal 

Transverse 

21 1 

jHRHH 
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204 


175 

.00130 

316 


173 


427 

119 

166 

.00141 

538 

114 

162 

.00148 

649 

109 

156 

.00149 

760 

103 

149 

.00156 

671 

97 

142 

.00160 


TABLE II. - DS MAR-M 246 STRESS- 
STRAIN PROPERTIES 
(LONGITUOINAL) 


plastic strain, 
percent 


0.1 

.2 

.4 

.6 

.6 

1.0 


Stress, MPa 

21 # C 


B9BI 

600 


875 

B30 


930 

850 


965 

855 

930 

970 

865 

945 

975 

870 

960 

980 
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Figure 1. - SSME high-pressure fuel 
turbopump first stage turbine blade. 
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Figure 2. - Turbine Inlet gas temperature for HPFTP. 
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Figure 5. - Distribution of heat transfer coefficients. 
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STRUCTURAL ANALYSIS DEMONSTRATION OF CONSTITUTIVE AND LIFE MODELS 


G. Bechtel, T.S. Cook, M.T. Tipton, and R.L. McKnlght 
General Electric Company 
Cincinnati, Ohio 


The overall objective of this program Is to demonstrate the applicability 
of NASA-developed advanced constitutive and life damage models for calculating 
cyclic structural response and crack Initiation In selected components of reus- 
able space propulsion systems. The computer model resulting from this program 
will enable the user to produce an accurate life prediction of hot gas path, 
life-limiting components of propulsion systems such as the space shuttle main 
engine (SSME). Program features Include: 

• Integration of Prior Developed Programs as Basic Approach 

Previously developed computer models addressing constitutive modeling 
and life damage will be combined In an advanced finite element anal- 
ysis to generate a sophisticated baseline life prediction program. 

• Test Verification of Synthesized Model 

A material data base will be established for the constitutive and 
life models parametrically Involving temperature, strain range, 
strain rate, mean strain/stress, and dwell time. Structural problems 
will be run Involving cyclic thermomechanical loading and multlaxlal 
stress states. The structural assurance models will be verified for 
multlaxlallty, cumulative damage, and damage mechanism Interaction. 

• Model Utilization With SSME Components 

The verified computer program will be used to accomplish the life 
predictions of three SSME critical components as evidence of the 
model functionality. 

Three engine components will be chosen for analysis. Only one of these 
components has been selected to date - a blade from the high pressure turbopump 
presently used In the SSME and shown In figure 1. The analysis will be carried 
out using the properties of PWA 1480, the proposed material for the next gener- 
ation of blades. Other candidates for the program are a redesigned blade and 
a nonrotating component. An analysis of these components will be carried out 
by using properties of appropriate materials that are of a quite different 
class than the highly anisotropic, single crystal structure of the PWA 1480. 
This will allow the analysis to test the selected constitutive and life models 
for quite different microstructures and symmetry classes. 

The constitutive model selected for analysis of the anisotropic PWA 1480 
Is the Stouffer-Dame model. However, the Walker model will be examined to some 
extent. For the Isotropic materials, the model developed by Ramaswamy will be 
employed, although others may be evaluated depending on data requirements and 
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what Is available. Various damage models will also be Investigated, namely 
continuous damage, hysteresis energy, and stralnrange partitioning. 



PLATFORM 


Figure 1 . 
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MECHANICAL PROPERTIES OF HIGH GRADIENT AND HOT ISOSTATICALLY PRESSED PWA 1480 


L.G. Frltzemeler and G.D. Schnlttgrund 
Rockwell International 
Canoga Park, California 


Single crystal superalloy (specifically, alloy 1480) turbine blades are 
now In development for application In the high pressure turbopump turbine sec- 
tions of the SSME. The mechanical and environmental demands placed on these 
components are quite severe and are significantly different from the require- 
ments of aircraft gas turbine applications for which the alloy was developed. 

In particular, the extreme thermal transients on start and shutdown and the 
very high mean stress vibratory loads lead to stronger fatigue excitations than 
generally encountered In gas turbine engines. In addition, the short design 
life, about 7.5 hours, of the SSME precludes creep and stress rupture as pri- 
mary failure mechanisms. 

It Is recognized that the mechanical properties of single crystal super- 
alloys, especially with respect to life under cyclic loads, are strongly 
dependent on mlcrostructural defects. Fatigue cracks In single crystal super- 
alloys have been found to Initiate at such mlcrostructural anomalies as Inclu- 
sions, secondary grains and, most frequently. Internal porosity. Control or 
elimination of these defects can be expected to significantly Increase the life 
of single crystal components. The goal of this program Is to Investigate the 
potential benefits due to reduction or elimination of casting porosity through 
high thermal gradient casting and hot Isostatic pressing applied to alloy 1480. 
Additionally, Improvement In fatigue life due to alternate heat treatment Is 
being evaluated. 

Internal porosity Is an Inherent feature of single crystal superalloy 
castings. The bulk of this porosity Is formed during solidification and Is 
caused by shrinkage and restricted fluid flow Into the Interstices of the den- 
drite arms. Increased casting thermal gradient can refine the microstructure 
through reduced dendrite arm spacing (DAS) and an attendant decrease In the 
size and density of the casting porosity. Internal casting porosity can be 
virtually eliminated by the proper application of hot Isostatic pressing (HIP) 
to the casting. The successful application of HIP to single crystal parts 
Involves overcoming a unique set of obstacles. Temperature control In produc- 
tion HIP vessels Is generally not tight enough to control within the as-cast 
solution heat treatment range of alloy 1480. Cooling rates from temperature 
are also not rapid enough to provide optimum properties. Post-HIP solution 
heat treatment cause pore resurgence through Klrkendall diffusion. Improper 
temperature and pressure application In the HIP cycle can also cause recrystal- 
lization around the closing pores. Very small amounts of free carbon In the 
HIP vessel cause surface carburization. With the successful removal of casting 
porosity, Improvements In fatigue life and Increased resistance to hydrogen 
environment embrittlement ( HEE ) may be obtainable through tailored heat treat- 
ments. Increased slip dispersal through optimum gamma prime distribution will 
provide Improvements In both fatigue and HEE behavior. 
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This program has been designed to evaluate the potential benefits to 
alloy 1480 material properties due to high thermal gradient casting, hot Iso- 
static pressing, and alternative heat treatment. Alloy 1480 castings have been 
obtained from vendors representing the extremes of commercial casting thermal 
gradients. Quantitative characterization of the DAS and pore distributions has 
been conducted. A HIP schedule, which avoids the many pitfalls of single crys- 
tal HIP has been devised for alloy 1480. In addition, an alternate heat treat- 
ment for alloy 1480 has been devised which provides benefits In mechanical 
properties. The standard gradient cast test material Is being evaluated In the 
standard heat treated condition as a baseline and In the HIP plus alternate 
heat treated condition. The high gradient cast material Is being evaluated In 
alternate heat treat and HIP plus alternate heat treated conditions. Evalua- 
tions focus on demonstration of the benefits due to the applied processes, 
especially In the area of cyclic life. This Interim period report will present 
current results and procedures derived from the program. 
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EFFECT OF HYDROGEN ON DEFORMATION STRUCTURE AND 
PROPERTIES OF CMSX-2 NICKEL-BASE SINGLE-CRYSTAL SUPERALLOY 

M. Dollar, I.M. Bernstein, S. Walston, 

F. Prinz, and A. Domnanovich 
Carnegie Mellon University 
Pittsburgh, Pennsylvania 

Material used in this study was a heat of the alloy CMSX-2 whose chemical 
composition is given in Table I. This nickel-based superalloy was provided in the 
form of [001] oriented single crystals, solutionized for 3 hrs at 1375°C. It was then 
usually heat treated as follows: 1050° C/16h/air cool 4 £50°C/48h/air cool. The 

resulting microstructure is dominated by cuboidal, ordered y' precipitates with a 
volume fraction of about 75% and an average size of 0.5 ^m (Figure 1). 

In previous studies (1,2) of hydrogen effects in the alloy we have shown that a 
high hydrogen fugacity environment can degrade room temperature mechanical 
properties. Hydrogen was introduced by cathodic potentiostatic charging in a mixture 
of two molten salts, 57% sodium bisulfate and 43% potassium bisulfate at 150 ° C for 
5 hrs. The most significant effects of such charging were on tensile behavior, where 
a 34% loss of homogeneous elongation was found and on push-pull fatigue behavior 
where a several-fold decrease of the accumulated plastic strain to fracture resulted, 
over a wide plastic strain range down to 0.02%. A detailed microscopic 
interpretation of these hydrogen-induced effects was limited by the fact that the 
hydrogen introduced by the molten salt charging resided primarily in a surface layer 
about 50 /im in depth. 

To better understand and control the role of hydrogen, high pressure, high 
temperature hydrogen charging was carried out at Sandia National Laboratory- 
Livermore, CA. This entailed gas phase charging for 15 days at 300 ° C under a 
pressure of 138MPa, and furnace cooling to room temperature in 16 hrs. under 
pressure. A resultant, uniform hydrogen concentration of the order of 5000 appm 
across the entire cross-section of 3 mm thick samples was obtained by such a 
procedure. 

As shown in Tables II and III, and in Figure 2, this charging mode more 
dramatically changes tensile and fatigue behavior than molten salt charging. The 
presence of this significant amount of dissolved hydrogen, while not significantly 
changing the relatively high yield strength, does promote significant work hardening, 
in contrast to hydrogen-free material where softening is found, and also dramatically 
reduces the true uniform elongation to almost zero. In fatigue, the number of cycles 
to failure now drops two order of magnitude and the occurrence of strain localization 
is more pronounced then in either the hydrogen-free case, or when the hydrogen was 
localized at the surface. 
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TEM and SEM evidence provides a picture of how hydrogen can affect the 
deformation structure and the fracture of CMSX-2 and can be used to correlate 
changes in mechanical properties with the presence of hydrogen in both the y and y 
phases. In brief, the most compelling hydrogen induced-changes in deformation 
structure are : i) enhanced dislocation accumulation in the y matrix; and ii) more 
extensive cross-slip of superdislocations in the y' precipitates. These are shown in 
Figures 3 and 4, respectively (for comparison see also Figures 5 and 6 which 
illustrate the deformation structure of H-free material). The enhanced dislocation 
density in y acts to decrease the mean free path of a superdislocation, while easier 
cross slip hinders superdislocation movement by providing pinning points in the form 
of sessile jobs. Both processes contribute to the increase of flow stress and the 
notable work hardening that occurs prior to fracture. 

Previously (1,2), examination of fractured tensile samples revealed that 
uncharged samples fractured predominantly along planes parallel to multiple, 
crystallographic {111} slip systems. The fractured surface was jagged and significant 
localized and general ductility was observed. In the present study, the hydrogen- 
charged samples were found to be much more brittle in appearance with a fracture 
surface that was quite planar and oriented normal to the tensile axis. A most 
striking SEM observation was that microcracks appear to follow <001> oriented y/y' 
interfaces (Figure 7) suggesting that preferrential cracking parallel to {001} planes 
occurs. However, TEM observations strongly suggests that a difference exists in the 
macroscopic and microscopic nature of hydrogen-induced crack initiation and 
propagation. In particular. Figure 8 clearly shows that cracking is associated with 

processes in the y, manifested as localized tearing in thin foils. While at first sight 

these appear to follow {001} planes, in reality they are comprised of short segments 
on multiple {111} planes. While tearing in thin foils is clearly not identical with 

microcracking in bulk material, it is suggestive and we believe accurately so, that 

embrittlement results from hydrogen-enhanced strain localization, resulting in 
dislocation intensification in the narrow ribbons of y and premature plastic failure 
due to strain exhaustion and the apparent mode I fracture. 

This work is being supported by NASA-Lewis under the technical direction of 
Dr. Robert Dreshfield. 
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TABLE I 


CMSX-2 Composition (wt. %) 


Cr 

8.0 

W 

8.0 

Ta 

6.0 

Al 

5.6 

Co 

4.6 

Ti 

1.07 

Mo 

0.6 

Fe 

0.08 

Si 

0.015 

Ni 

balance 

C 

15 ppm 

S 

10 ppm 

N 

4 ppm 

0 

2 ppm 


TABLE II 

Tensile Behavior as a Function of Charging Condition 


Charging 

Conditions 

YS (MPa) 
(0.2% offset) 

UTS (MPa) 

Homogeneous 
Elongation (%) 

Elongation 3 

Loss (%) | 

uncharged 

917 

987 

21.7 

1 

1 

1 

molten salt 
5h at 150°C 

856 

976 

14.4 

33.6 

gas phase 

5 days at 
1 38MPa & 300 °C 

935 

1123 

1.2 

94.5 


111 






















TABLE III 


Fatigue Behavior as a Function of Plastic Strain Range 
and Charging Conditions 







Accum. Plastic Strain | 

A€p(%) 

Charging 

Conditions 

Cycles to Failure 
uncharged charged 

n u /n c 

(n x 
uncharged 

AC P ) (%> 
charged 

0.4 


57 

16 

3.6 

46 

13 

0.2 

molten salt 

131 

57 

2.3 

52 

23 

0.1 

5h at 150°C 

214 

80 

2.7 

43 

16 

0.02 


3348 

1260 

2.7 

137 

53 

0.02 


6664 

21 

317 

267 

.84 

0.02 

gas phase 

3085 

16 

193 

123 

.56 

0.02 

15 days at 
138MPa & 300 °C 

2483 

29 

86 

99 

1.16 


R = -1. Room Temperature Testing 


112 









Engineering Stress (MPa) 



-a- Uncharged 
Molten Salt 
ho- Gas Phase 










c 


eformation structure 
SX-2; e = 4 percent. 


ure Figure 6. - Typical deformation structure 
of hydrogen- free CMSX-2; c = 5 percent; 
(foil cut parallel to (111) slip plane). 


ORIGINAL PAGE IS 
OF POOR QUALITY 


114 









ORIGINAL PAGE IS 
OF POOR QUALITY 


Figure 7. - Fracture surface of hydrogen- charged CMSX-2. 


g in a 
was hy 


in foil cut close to fracture sur- 
rogen charged . ) 

115 


N87 - 22783 


PROBABILISTIC STRUCTURAL ANALYSIS TO EVALUATE THE STRUCTURAL DURABILITY OF 

SSME CRITICAL COMPONENTS 

Christos C. Chamis 
NASA Lewis Research Center 
Cleveland, Ohio 


It is becoming increasingly evident that deterministic structural analysis 
methods will not be sufficient to design critical structural components for 
upgraded space shuttle main engines (SSME). Structural components in the SSME 
are subjected to a variety of complex (cyclic and transient) loading conditions 
including high temperatures and severe thermal gradients. Most of these are 
quantifiable only as best engineering estimates. These complex loading condi- 
tions subject the component material to coupled nonlinear behavior which 
depends on stress, temperature, and time. The nonuniform nature of this 
coupled nonlinear material behavior makes it very difficult to characterize 
experimentally and perhaps impossible to describe deterministically. 

In addition, critical SSME structural components are relatively small. 
Fabrication tolerances on these components, which in essence are small thick- 
ness variations, can have significant effects on the component structural 
response. Fabrication tolerances by their very nature are statistical. 
Furthermore, the attachment of the components to the structural system gener- 
ally differs by some indeterminant degree from what is assumed for design pur- 
poses. In summary, all four fundamental aspects - (1) loading conditions, (2) 
material behavior, (3) geometric configuration, and (4) attachment - on which 
structural analyses are based, are of statistical nature. One direct way to 
formally account for all these statistical aspects is to develop probabilistic 
structural analysis methods where all participating variables are described by 
appropriate probabilistic functions. 

NASA Lewis Research Center is currently developing probabilistic struc- 
tural analysis methods for select SSME structural components. Briefly, the 
deterministic, three-dimensional, inelastic analysis methodology developed 
under the Hot Section Technology ((HOST) and R&T Base Programs) is being aug- 
mented to accommodate the complex probabilistic loading spectra, the thermo- 
viscoplastic material behavior, and the material degradation associated with 
the environment of space propulsion system structural components representative 
of the SSME such as turbine blades, transfer ducts, and liquid-oxygen posts 

(fig. 1). 

The development of probabilistic structural analysis methodology consists 
of the following program elements: (1) composite load spectra, (2) probabil- 

istic structural analysis methods, (3) probabilistic finite element theory - 
new variational principles, and (4) probabilistic structural analysis applica- 
tion. In addition, the program includes deterministic analysis elements: (1) 

development of structural tailoring computer codes (SSME/STAEBL) , (2) develop- 
ment of dynamic creep buckling/ratcheting theory, (3) evaluation of the dynamic 
characteristics of single-crystal SSME blades, (4) development of SSME blade 
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damper technology, and (5) development of integrated boundary elements for hot- 
fluid structure interaction. 

The SSME Structural Durability Program is multidisciplinary integrated and 
is a joint effort of NASA Lewis in-house research, contracts, grants, and sup- 
port service personnel (table I). The research activities in the various pro- 
gram elements have collectively led to significant technical progress. The 
objective of the presentations included in this session is to summarize the 
highlights of this technical progress. 


TABLE I. - COMPUTER CODES DEVELOPMENT SUMMARY 


Code 

Description 

Status 

NESSUS 

Nonlinear evaluations of stochastic 
structures under stress 

Initial version 
operational 

CLS 

Composite load spectra probabilistic 
load simulation 

Initial version 
operational 

SSME/STAEBL 

Structural tailoring of SSME turbopump 
blade 

Code available for 
release 

VATFEL 

Variational theory for probabilistic 
finite elements 

Near completion 

DYCREBURT 

Dynamic creep buckling/ratcheting theory 

Near completion 

BEFSIN 

Boundary elements for fluid structure 
interaction 

Planning phase 

NL-COBSTRAN 

Nonlinear structural analysis for high- 
temperature metal matrix composite turbine 
blades 

Operational 
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Figure 1 . - PSAM will be initially developed for select SSME structural components. 
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PROBABILISTIC STRUCTURAL ANALYSIS METHODS FOR SELECT SPACE 


PROPULSION SYSTEM STRUCTURAL COMPONENTS (PSAM) 


T.A. Cruse 

Southwest Research Institute 
San Antonio, Texas 


The objective of this 5-year contract effort 1 Is the development of 
several modular structural analysis packages capable of predicting the proba- 
bilistic response distribution for key structural variables such as maximum 
stress, natural frequencies, transient response, etc. The structural analysis 
packages are to Include stochastic modeling of loads, material properties, 
geometry (tolerances), and boundary conditions. As shown In figure 1, the 
solution Is to be In terms of the cumulative probability of exceedance distri- 
bution (CDF) and confidence bounds. Two methods of probability modeling are 
to be Included as well as three types of structural models - probabilistic 
finite-element method (PFEM); probabilistic approximate analysis methods 
(PAAM); and probabilistic boundary element methods ( PBEM) . 

The purpose In doing probabilistic structural analysis Is to provide the 
designer with a more realistic ability to assess the Importance of uncertainty 
In the response of a high performance structure. PSAM tools will estimate 
structural safety and reliability, while providing the engineer with Informa 
tlon on the confidence that should be given to the predicted behavior. Perhaps 
most critically, the PSAM results will directly provide Information on the 
sensitivity of the design response to those variables which are seen to be 
uncertain. Design situations where the performance of the structure strongly 
depends on variables whose values are uncertain can be Identified and correc- 
tive steps taken. 

The PSAM tools will reflect three levels of user-defined uncertainties, 

as shown In figure 2. The first level Is one In which the uncertainty Is from 

part to part or application to application. Thus, for example, one turbine 
blade has a different fundamental frequency from a second, randomly selected, 
blade. The second level of uncertainty reflects the fact that the statistics 
on loading or material behavior or tolerances may vary from one region of the 
structure to another. The third, and most detailed, level of uncertainty Is 
within the local behavior of a continuum. In this level, for example, the 
strains within a structure are to be randomly correlated to the displacement 
derivatives; the stresses are to be randomly correlated with the strains; and 
so on. 

Three major analysis codes are to be developed following the schedule 
shown In figure 3. The first of these, the PFEM code, has been under develop- 
ment for just over 2 years. A linear version of that code has been developed 

(fig. 4). Validation studies have been Initiated and the code will be applied 
this year to an SSME turbine blade analysis. The PAAM modeling effort will not 


^ork done under NASA Contract NAS3-24389. 


{’RECEDES 




v i:0 


121 


be Initiated until FY 88. Following a 1-year study of advanced analysis 
methods, the development of a probabilistic boundary element code has been 
Initiated. The PBEM code structure will follow that being established for the 
PFEM code NESSUS. 

The current NESSUS (Numerical Evaluation of Stochastic Structures Under 
Stress) code structure Is shown In figure 4. The code currently consists of 
four major modules: NESSUS/EXPERT, NESSUS/FPI, NESSUS/PRE, and NESSUS/FEM. 

The EXPERT module provides a user Interface to NESSUS to facilitate the 
development and Interpretation of probabilistic data. Fast Probability Inte- 
gration (FPI) Is the currently used reliability estimation algorithm for 
assessing response distributions and confidence levels. The PRE module reduces 
random fields with partial correlations to a set of Independent random field 
variables using modal analysis methods. The finite-element method (FEM) module 
provides a rapid solution method especially developed to perform sensitivity 
analysis of a structural model. Two subsequent presentations In this same 
volume will discuss the FPI and EXPERT modules, and then the PRE and FEM 
modules . 

Figure 5 Illustrates the team approach being used In the PSAM contract 
effort. SwRI Is responsible for the system Integration and code validation 
efforts, as well as the development of the FPI and EXPERT modules. MARC Ana- 
lysis Research Corporation Is responsible for the FEM and PRE modules. Rocket- 
dyne Division has defined four verification problems and will apply the NESSUS 
codes to the probabilistic modeling of these applications. Their presentation 
follows those of SwRI and MARC. Professors Wlrschlng (University of Arizona) 
and Atlurl (Georgia Institute of Technology) and their graduate students are 
providing support on the development of advanced reliability and finite-element 
methods, respectively. Their presentations conclude the set of PSAM 
presentations . 

The goal of the PSAM contract Is being attained. A clear methodology for 
the probabilistic simulation of structural response Is established. A linear 
modeling capability for static and vibrating systems exists and has been for- 
warded to NASA. The current effort Is focusing on a fuller Integration of 
NESSUS/FEM and NESSUS/FPI algorithms In order to more accurately and effi- 
ciently define the tails of the distributions. MARC Is developing nonlinear 
solution procedures for probabilistic structures, as well as Improving the FEM 
technology being used. The resulting structural analysis methods are seen to 
be capable of providing new and very significant Insight Into the reliable 
performance of complex structural systems. 
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PSAM INTEGRATES FEM TECHNOLOGY 
WITH PROBABILISTIC METHODS 
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THE NESSUS FINITE- ELEMENT CODE* 

J.B. Dias, J.C. Nagtegaal, and S. Nakazawa 
MARC Analysis Research Corporation 
Palo Alto, California 


The NESSUS finite element code is being developed by MARC Analysis Research 
Corporation as part of the probabilistic structural analysis (PSAM) effort, 
coordinated by Southwest Research Institute for the NASA Lewis Research Center. 
The objective of this development is to provide a new analysis tool which 
integrates the structural modeling versatility of a modern finite element code 
with the latest advances in the area of probabilistic modeling and structural 
reliability. Version 2.0 of the NESSUS finite element code was released to the 
members of the PSAM team last February, and is currently being exercised on a set 
of problems which are representative of typical SSME applications. NESSUS 2.0 
allows linear elastostatic and eigenvalue analysis of structures with uncertain 
geometry, material properties and boundary conditions, which are subjected to a 
random mechanical and thermal loading environment. 

The NESSUS finite element code is a key component in a broader software 
system consisting of five major modules. NESSUS/EXPERT is an expert system under 
development at Southwest Research Institute, with the objective of centralizing 
all component-specific knowledge useful for conducting probabilistic analysis of 
typical SSME components. NESSUS/FEM contains the finite element code used for 
the structural analysis and parameter sensitivity evaluation of these components. 
The task of parametrizing a finite element mesh in terms of the random variables 
present is facilitated with the use of the probabilistic data preprocessor in 
NESSUS/PRE. An external database file is used for managing the bulk of the data 
generated by NESSUS/FEM. To complete the analysis, Southwest Research Institute 
has developed the probabilistic analysis module NESSUS/FPI, which extracts from 
the database the information needed for generating probability distributions for 
the desired response variables. 

Probabilistic finite element analysis involves the computation of the 
effects of small variations of the random variables on the overall response of 
the structure. Since probabilistic models of realistic structures often require 
large finite element models parameterized by many random variables, this step is 
very computation-intensive and will tend to govern the solution cost. The 
strategy adopted in the NESSUS finite element code is based on an iterative 
algorithm akin to the modified Newton method. With this approach, all 
computations on the perturbed system can be performed at the element level, 
resulting in substantial savings on memory requirements. This approach also 
eliminates the need to compute and store explicit partial derivatives of the 
element matrices and vectors, or to assemble these into additional global arrays. 
As a result, the perturbation of the element parameters can be made independent 
of element formulation, allowing the introduction of the newest element 
technology into the existing algorithm with a minimal amount of new software 
development. Significant improvements in convergence speed have been 


*Work performed under NASA contract NAS3- 24389. 
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demonstrated with the use of modern algorithms for nonlinear programming, which 
have been implemented in NESSUS within the past year. 

The iterative algorithms implemented in the NESSUS finite element code can 
be shown to be inherently consistent, but must satisfy known stability conditions 
in order to achieve convergence. In most cases this will not present a problem, 
since loss of stability will only occur for very large stiffness changes, which 
are well beyond the range of uncertainty observed in the physical structure. 
However, there are ill-conditioned problems for which even seemingly small 
perturbations will engender large changes in some stiffness terms. This is often 
the case when constraint equations are present in the problem formulation. Such 
problems are typically encountered in the analysis of strongly anisotropic or 
incompressible materials (e.g. composites, rubber mounts and seals), deviatoric 
rate-independent plasticity, or the thin limit of the Reissner-Mindlin shell 
theory. A solution to this problem will involve the development of smart 
adaptive algorithms for selecting the perturbation size in order to ensure 
convergence of the algorithm. 

Significant progress was made over the last year on issues related to the 
modeling of probabilistic structures and the development of a good user 
interface. It was soon realized that classical shell theories lack the 
generality needed to capture localized effects which are often important in the 
analysis of SSME-related components. This resulted in the development of a new 
family of finite elements based on an assumed strain formulation and designed to 
achieve superior performance in bending problems. These elements allow the 
definition of continuous pressure fields at the surface nodes using a new 
boundary identification algorithm. Current development efforts are aimed at 
enhancing the element's thermal strain response and extending the new formulation 
to problems involving material nonlinearity. 

The external database file is used for retaining a permanent record of the 
parameter sensitivity data generated by the NESSUS code. This database resides 
in a binary (unformatted) direct-access file, and is structured as a two-way 
ordered linked list. This type of data structure allows the insertion, deletion 
and replacement of individual entries without the need to move large blocks of 
data. It is therefore possible to maintain and expand an existing database with 
additional results obtained in the course of multiple runs of the NESSUS finite 
element code. Each database file is self-contained, allowing direct access to 
the perturbation data from within other modules in the NESSUS system. 

Planned developments for the coming year include extending the capabilities 
of the NESSUS code to perform probabilistic analysis in situations involving 
material nonlinearity and transient dynamics. This will require tracking 
separate time-histories for all perturbed problems and transferring perturbation 
data between increments. The ability to conduct random vibration analysis on 
uncertain structures is also planned for the near future. The development of 
probabilistic finite deformation analysis in NESSUS is currently planned for the 
fourth year of the PSAM effort. 
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FIGURE Is AN OVERVIEW OF THE NESSUS CODE 




FIGURE 2: UNCOUPLED RANDOM VARIABLES USED TO CHARACTERIZE GEOMETRICAL 

UNCERTAINTY MAY BE IDENTIFIED WITH SIMPLE GEOMETRY IMPERFECTIONS 
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FIGURE 


Unperturbed Perturbed 



FIGURE 3: GRAPHICAL INTERPRETATION OF MODIFIED NEWTON ITERATION 


Unperturbed 



: STABILITY BOUNDS DEFINE THE ALLOWABLE RANGE OF PERTURBATION SIZES 

WHICH ARE CONVERGENT WITH THE MODIFIED NEWTON ALGORITHM 
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FIGURE 5: LOCALIZED THERMAL AND MECHANICAL EFFECTS CRUCIAL TO THE ANALYSIS OF 

COMPONENTS SUCH AS THIS BLISTER SPECIMEN MODEL CANNOT BE ADEQUATELY CAPTURED 

USING SIMPLIFIED SHELL-BASED MODELS 



Increment 0 
Perturbation 0 

Increment 0 
Perturbation 1 

Increment 0 
Perturbation 2 

Increment 1 
Perturbation 0 

Increment 1 
Perturbation 1 

Increment 1 
Perturbation 2 

Increment 3 
Perturbation 0 


FIGURE 6: AN EXTERNAL DATABASE FILE IS USED TO RETAIN A PERMANENT RECORD OF 

THE PARAMETER SENSITIVITY DATA GENERATED BY THE NESSUS CODE 
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PROBABILISTIC STRUCTURAL ANALYSIS VERIFICATION STUDIES 


K.R. Rajagopal 

Rocketdyne Division, Rockwell International 
Canoga Park, California 


The space propulsion system components are subjected to complex mechanical 
and thermal loading. Typically, these components are analyzed using the 
state of the art deterministic methodologies in linear and nonlinear 
structural analysis domain. The components are then designed using a 
factor of safety approach. The purpose for doing a probabilistic 
structural analysis of space propulsion system components has been 
explained in one of the earlier papers on PSAM appearing in this same 
volume. 

During the first year and part of the second year of the contract, 
identification and definition of environments and analysis techniques that 
are used to analyze the typical space propulsion system components were 
discussed and presented (Ref. 1). The basic objective of the verification 
efforts is to apply the probabilistic structural analysis methods 
developed and implemented in NESSUS (Numerical Evaluation of Stochastic 
Structures Under Stress) code to typical space propulsion components 
(Figure 1). The chosen typical components are turbine blade, high 
pressure duct, Lox post, and transfer tube liners. Since analysis options 
of increasing levels of sophistication are implemented in NESSUS 
incrementally, the verification efforts are also tailored to have 
increasing levels of sophistication during the progression of the 
contract. The current released version of the code is limited to linear 
structural analysis. 

Considering turbine blade component first, several simple verification 
studies were conducted. The studies exercised the solid element type, 
typical random variables, and various solution strategies. The results of 
the simple verification studies (Ref. 2) aided in establishing confidence 
in the code, identified its limitations in user interface, finite element 
and analysis deficiencies. The studies conducted on simple models so far 
indicate that for the chosen random variables and their normal range of 
perturbation found in practice, the implemented solution strategies were 
adequate. On the other hand, the studies also showed that when the 
perturbations were large, in the context of doing a sensitivity analysis 
over a wide range, the implemented solution strategies did not converge, 
and in some cases divergence occurred. The divergence or slow convergence 
issue has been addressed in the development of fully integrated 
probabilistic structural analysis package that causes the reformulation of 
the system equations as needed for convergence. 

The studies included geometry perturbation (variations in tilt angle), 
material orientation perturbation (+10°), material property perturbation 
and load (centrifugal load, pressure, temperature) perturbations. The 
studies resulted in improved user interfaces and hastened the 
implementation of more advanced solid elements that behave better in a 
variety of stress fields including temperature gradients. 
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Several simple verification studies have been conducted on frequency 
extraction of the perturbed structure using perturbation technique, which 
is of special interest in turbine blade components. Based on the review 
of results, strategies for improvements in accuracy and a reduction in 
computing time to obtain the frequencies of perturbed structure are being 
sought. 

A solid finite element turbine blade model, representative of blades found 
in space propulsion system components, has been generated and is currently 
being exercised (Figure 3). A basic set of random variables, consistent 
with first cut design requirements, have been identified (Figure 3). 
These random variables will be used in the study to determine the 
probability density function (PDF) of stress at critical areas and the 
PDF of frequencies. Some of the other important random variables such as 
variables contributing to the damper effectiveness, dynamic pressure, 
support conditions and transient conditions that are necessary for the 
final design are not considered in the initial exercise. 

Following the turbine blade analysis, verification efforts on other 
typical components such as high pressure ducts (Figure 4), Lox post 
(Figure 5), and transfer tube liners (Figure 6) will follow. The random 
variables and solution types have been carefully chosen to exercise and 
verify a wide range of analysis options and element types that are being 
developed in NESSUS and also in determining the probabilistic response of 
primary component specific design parameters. 
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NESSUS/EXPERT AND NESSUS/FPI IN THE PROBABILISTIC STRUCTURAL 
ANALYSIS METHODS (PSAM) PROGRAM 

O.H. Burnside 

Southwest Research Institute 
San Antonio, Texas 

NESSUS (Numerical Evaluation of Stochastic Structures under Stress) is 
the primary computer code being developed in the NASA PSAM project. It 
consists of four modules NESSUS/EXPERT, NESSUS/FPI, NESSUS/PRE and 
NESSUS/FEM. This presentation concentrates on EXPERT and FPI, while PRE and 
FEM are discussed in another presentation. 

One challenge of the PSAM effort is the effective integration of advanced 
finite element and probabilistic methods. A code with linear static and 
dynamic capabilities has been provided to NASA. However, it is clear that, in 
the final version, the user must be provided with an interface program to 
effectively use features such as nonlinear analysis and confidence band 
estimation. Such an interface program will also expedite the process of 
conducting the analyses necessary for NESSUS verification. 

To provide an effective interface between NESSUS and the user, an expert 
system module called NESSUS/EXPERT is being developed. That system uses the 
CLIPS artificial intelligence code developed at NASA-JSC. The code is 
compatible with FORTRAN, the standard language for codes in PSAM. The user 
interacts with the CLIPS inference engine, which is linked to the knowledge 
database as shown in Figure 1 . 

The essential features of EXPERT are its automated user input and 
automated results. The EXPERT module will provide the user with features such 
as interactive HELP, data set consistency checking, and defaults for the 
statistical models of random variables. EXPERT will also assist the analyst 
in managing the large database produced in perturbing the random variables. 
Such variables may include material properties, geometry, boundary conditions 
and loading. Because of the potentially large number of random variables, 
this process must be automated to free the analyst from basically a 
bookkeeping task. For analysis of certain critical SSME components, EXPERT 
will choose component specific perturbations. For example, in the case of a 
turbine blade shown in Figure 2, perturbations may be applied directly to 
parameters defining the blades geometry such as the twist and tilt angles and 
blade thickness. 

The perturbation database generated by NESSUS/FEM and managed in EXPERT 
is used to develop the so-called response or performance model in the random 
variables. Figure 3 illustrates such a model in which natural frequency is a 
random function of material modulus. It is from this performance model that 
the probabilistic response is computed. Two independent probabilistic methods 
are available in PSAM for the computation of the probabilistic structural 
response. These are the Fast Probability Integration (FPI) method and Monte 
Carlo simulation. FPI is classified as an advanced reliability method and has 
been developed over the past ten years by researchers addressing the 
reliability of civil engineering structures. Monte Carlo is a well- 
established technique for computing probabilities by conducting a number of 
deterministic analyses with specified input distributional information. 
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For structural systems, the probability of failure is generally low. In 
such situations, Monte Carlo is inefficient since a large number of 
simulations are required to confidently predict low probability levels. The 
efficiency of FPI , on the other hand, is not tied to the probability level and 
can accurately predict the tails of the response distribution. However, Monte 
Carlo is still useful in PSAM to check the accuracy and robustness of the FPI 
algorithms. For a given performance function, FPI can compute point 
probability estimates or obtain full distributional information in terms of 
the cumulative probability distribution. In general, the performance function 
may be nonlinear in the random variables and contain mixed probability 
distribution types. Figure 4 qualitatively compares typical FPI and Monte 
Carlo results and illustrates the increasing inaccuracy in Monte Carlo at the 
low probability levels. 

The process used by FPI to make the probability estimates may be 
considered a problem of constrained minimization. This is illustrated in 
Figure 5. Let us assume that material density and modulus of elasticity are 
random variables in an eigen-frequency analysis. The response function at a 
particular frequency is shown, with the lines of constant probability given by 
the circles. The design or most probable point is defined as the point with 
the highest joint distribution at a given value of the response function. The 
design point is geometrically located at the minimum distance 8, called the 
safety index. In a first-order reliability analysis, 8 is related to the 
probability of exceeding the specified frequency. 

The FPI algorithm uses an efficient iteration technique to converge to 
the most probable point. The first estimate of the design points uses the 
NESSUS-generated database at the mean state. From this initial guess, the 
automated algorithms in FPI use successive NESSUS perturbations at other 
points to converge to the design point as illustrated in Figure 6. 
Convergence is usually obtained in several (three to five) iterations. This 
procedure gives the probability of exceedence at a particular value of the 
response function. This is called a point probability estimate. Cumulative 
distributional information can be obtained by running FPI at several values of 
the response variables as shown on Figure 4. 

In summary, the NESSUS/EXPERT and NESSUS/FPI modules provide the 
designer/analyst with computational methodologies and software for effectively 
and efficiently evaluating design sensitivities and uncertainties and 
quantifying the structural performance of SSME components. 
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EXPERT KNOWLEDGE WILL GOVERN 
USER INTERFACE 



FIGURE 1. 

NESSUS PERTURBATIONS TO BE AUTOMATED 



FIGURE 2. 


MATERIAL PROPERTIES 

• ANISOTROPIC CONSTANTS 

• ANISOTROPIC ORIENTATIONS 

• STRESS-STRAIN CURVES 

GEOMETRY 

• TWIST 

• TILT 

• THICKNESS 

• ETC. 

BOUNDARY CONDITIONS 

• LOADS 

• CONSTRAINTS 


141 




FPI CODE RELIES ON RESPONSE MODEL 


INPUT DATA 


HI SI’ONSI MODI I 


CDF 



+ FREQUENCY 


MODULUS 



• 100 RANDOM VARIABLES 

• INDEPENDENT STATISTICS 

- WEIBULL 

- LOG NORMAL 

- NORMAL 

- EMPIRICAL 

FREQUENCY 

FIGURE 3. 

FPI CODE VALIDATED BY MONTE CARLO 



FPI ATTRIBUTES 

• UP TO 100 TIMES FASTER THAN MC 

• EFFICIENCY INDEPENDENT OF PROBABILITY 
LEVEL 

• CAN USE LINEAR/QUADRATIC DESIGN 
SENSITIVITY DATA 

• CAN BE COUPLED TO FEM 

• USES ANY INPUT STATISTICS 

• COMPUTES RISK LEVEL 

• ESTIMATES CONFIDENCE BANDS 
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FPI ALGORITHM APPROXIMATES 
NONNORMAL VARIABLES 



• APPROXIMATIONS BASF3D ON 
EQUAL CDF NFAR DFSIGN 
POINT 

• AREA OF EXCEEDANCE 
CALCULATED FOR EQUIVALENT 
NORMAL DISTRIBUTIONS 

• ERRORS GREATEST AWAY 
FROM DESIGN POINT 

• ERROR PROBABILITIES ARE 
SMALL 

• PROBABILITY GIVEN FROM 
(3 DISTANCE 


FIGURE 5. 


FPI ALGORITHM ITERATES TO MAXIMUM 
LIKELIHOOD POINT 



• ESTIMATE LINEAR RESPONSE 
FUNCTION 

• CALCULATE MOST LIKELY 
DESIGN POINT 

• SHIFT RESPONSE 
FUNCTION SURFACE 

• RECOMPUTE DESIGN 
POINT 


FIGURE 6. 
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ADVANCED PROBABILISTIC METHOD DEVELOPMENT 

P.H. Wlrschlng 
The University of Arizona 
Tucson, Arizona 


Advanced structural reliability methods are utilized on the PSAM project 
to provide a tool for analysis and design of space propulsion system 
hardware. The role of the effort at the University of Arizona is to provide 
reliability technology support to this project. Reliability refers, in 
general, to application of probabilistic and statistical methods to manage 
uncertainties in engineering design. 

PSAM computer programs will provide a design tool for analyzing uncer- 
tainty associated with thermal and mechanical loading, material behavior, 
geometry, and the analysis methods used. Specifically, reliability methods 
are employed (a) to perform sensitivity analyses (a study of the importance of 
each design factor), (b) to establish the distribution of a critical response 
variable, (e.g., stress, deflection), (c) to perform reliability assessment, 
and (d) ultimately to produce a design which will minimize cost and/or weight 
(Fig. 1). 

Uncertainties in the design factors of space propulsion hardware, as 
defined in Fig. 1, are described by probability models constructed using 
statistical analysis of data. Two examples are provided in Fig. 2 and Fig. 

3. Turbine blade twist (a geometric variable) has a statistical distribution 
(Fig. 2), and significant variability in material properties are observed in 
laboratory tests (Fig. 3). Statistical methods are employed to produce a 
probability model, i.e., a statistical synthesis or summary of each design 
variable (e.g., each in Fig. 4) in a format suitable for reliability 
analysis and ultimately, design decisions. 

Given the probability model for each random variable, and the design 
equation, i.e., how the variables are related (e.g., through the NESSUS/FEM 
computer code), two basic methods for computing the distribution of a response 
variable are available. They are, fast probability integration (FPI ) , (Fig. 

4) and Monte Carlo (Fig. 5). FPI is a complex numerical algorithm which 
provides "fast" probability calculations, accurate to about ± 5% (See Refs. 1 
and 2). 

The NESSUS generated perturbed data set, an example of which is shown in 
Fig. 5, is used as a function of the input variables. Monte Carlo produces an 
exact solution, but it is expensive to run. There are a large number of 
different Monte Carlo schemes available; methods of improving efficiency of 
Monte Carlo, tailored for PSAM application, are being studied. Availability 
of the response function shown in Fig. 5 is a requirement for performing a 
fast Monte Carlo analysis. The NESSUS-generated, approximate response 
function can be used to replace the "exact" simulation models needed by Monte 
Carlo methods. 
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The result of PSAM analysis is a distribution of a response variable. But 
for design purposes, it is often important to be able to specify also the 
confidence associated with a probability statement. Confidence Intervals 
provide a quantitative description of the quality of the results of the 
probabilistic analyses (Fig. 6). Efficient methods for constructing confi- 
dence intervals are being studied. 

In summary, the University of Arizona effort on PSAM is to provide 
technical support in statistical characterization of design variables, 
advanced reliability analysis methods, and confidence interval construction 
(Fig. 7). 
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RELIABILITY METHODS ARE EFFECTIVE 
TOOLS FOR DESIGNERS 


SOURCES OF UNCERTAINTY: 

• LOADING 

• MATERIAL BEHAVIOR 

• GEOMETRY 

• STRUCTURAL MODELING 

• ANALYSIS METHOD 

• SMALL SAMPLE SIZE 



DESIGN APPLICATIONS 

• SENSITIVITY ANALYSIS 

• VARIABILITY ASSESSMENT 

• RELIABILITY ESTIMATION 

• COST/WEIGHT REDUCTION 

• CONFIDENCE LEVELS 


FIGURE 1. 
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STATISTICAL METHODS TRANSLATE 
DATA IN MODELS 



FIGURE 2. 


MATERIAL PROPERTIES EXHIBIT 
SIGNIFICANT VARIABILITY 



FIGURE 3. 
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FPI ALGORITHM IMPROVEMENTS 
DEMONSTRATED 


FATIGUE DAMAG E PR EDICTION MODEL 


• ENHANCED DISTRIBUTION LIBRARY 

- MATERIAL BEHAVIOR 

- LOADING UNCERTAINTY 

• IMPROVED ACCURACY 

• VALIDATION PROBLEMS 

- HIGHLY NONLINEAR 

- MANY VARIABLES 

• ACCURACY OF FPI i 5% 


[ X 2 1 X 2 
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EACH X; HAS DIFFERENT DISTRIBUTION 



• RELATIVE COMPUTER TIME 
FPI = 1 

MONTE CARLO = 60 
(CONVENTIONAL) 


FIGURE 4. 


FAST MONTE CARLO ALGORITHM (S) 
ESTABLISHED 


STRESS AT 
CRITICAL POINT 


• AN INDEPENDENT CHECK ON FPI 

• CONVERGES TO EXACT SOLUTION 



• MORE EXPENSIVE THAN FPI 

• USES NESSUS DATA BASE 


GEOMETRY 

VARIABLE 


• SEVERAL MONTE CARLO METHODS AVAILABLE 

• RELATIVE MERITS OF EACH METHOD UNDER 
STUDY 


MATERIAL 

BEHAVIOR 

VARIABLE 


NESSUS GENERATED DATA 

• DEFINES RESPONSE SURFACE 

• SUPPORTS MONTE CARLO SIMULATION 


FIGURE 5. 
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WHAT IS OUR CONFIDENCE IN THE 
ANSWER? 
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• DIFFICULT ANALYSIS PROBLEM 

• APPROXIMATE METHOD DEFINED 



• USES NESSUS DATA BASE 

• GIVES INTERVAL ESTIMATES 
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- VARIABILITY 


UA PROVIDES RELIABILITY 
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• STATISTICAL METHODS TO PRODUCE 
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- GEOMETRY 
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FIGURE 7. 
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STOCHASTIC AND HYBRID-STRESS PLATE/SHELL FINITE ELEMENTS FOR 

HOT-SECTION COMPONENTS 
S . N . At 1 u r i 

Georgia Institute of Technology 
Atlanta, Georgia 

The research effort in the Center for the Advancement of 
Computational Mechanics at Georgia Tech has two main thrusts. 
The first of these is the development of special approaches for 
the numerical stress analysis of solids and structures whose 
material and geometric properties are uncertain. The second 
seeks to develop and implement high-efficiency plate and shell 
e 1 emen t s . 

For a structure in which all loads, material properties and 
boundary conditions are certain, i.e. possess no random nature, 
it is possible to solve for precise values of displacement, 
stress and strain. On the other hand, for a structure in which 
material properties are uncertain, it is necessary to describe 
the structural response, e.g. displacement, in a manner which 
both incorporates and quantifies the resulting uncertainty. To 
this end the finite element methodology has been r e - f o r mu 1 a t e d so 
as to allow for the determination of a statistical description of 
stress, strain, etc. in terms of their means, variances and 
intercorrelations . 

The stochastic element method, currently being implemented, 
will be able to more accurately portray the probabilistic nature 
of stress, strain, and displacement in actual structures. 

In analyses of plate— and shell— like structures which 
exhibit strong gradients in stress, strain, or temperature 
through the thickness, accuracy demands that so-called 
’degenerate' 3-D elements be used, as opposed to conventional 
plate or shell elements. The degenerate plate/shell elements 

currently available, however, admit aspect ratios only up to 
about 10 to 1 before their performance begins to degrade. As a 
consequence, use of degenerate shell elements has heretofore been 
too costly, except for simplified problems. 

Current research has provided a hy b r i d - s t r e s s shell element 
whose behavior is acceptable for aspect ratios as high as 30 to 1. 
Thus, substantially more complex analyses will be practicable as 
soon as this element is fully implemented. 

An additional advantage of the hybrid approach is that it 
permits more accurate s t r e s s - r eco ve r y at the upper and lower 
surfaces of the shell, an important consideration in high 
thickness-gradient applications. 
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The software associated with the above research is being 
implemented in the form of extensions to the Nessus code, 
developed jointly by Southwest Research Institute and MARC 
Analysis Corporation. The implementation is being carried out in 
the Center for the Advancement of Computational Mechanics at 
Georgia Tech. The hybrid shell element has been successfully 
tested in several small-deformation elastic analyses; 
implementation of the remaining capabilities is underway. The 
theoretical formulation of the stochastic elements is essentially 
complete; its implementation is just beginning. 

A c k now 1 e dgemen t s : The contributions to this effort from Messrs. 

K~r L. Chen and T. Fitzgerald, and Drs. K. W. Reed and C.-T. Yang 
are sincerely acknowledged. 
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PROBABILISTIC FINITE ELEMENTS* 


Ted Belytschko and Wing Kam Liu 
Northwestern University 
Evanston , III ino is 


In the design of critical components, such as rocket engine parts, we 
often find large uncertainties in material properties and loads, particularly 
when considering the ultimate capacity of the components. One reason for this 
is that the extreme conditions, such as temperature or load, are not known, 
nor are the material properties at extreme temperatures as well established. 
Therefore, a single finite element analysis of a component may be quite 
misleading since it gives no information on the range of responses that can be 
expected. 

Although analysts often try to guard against this shortcoming by varying 
several of the parameters either arbitrarily or on the basis of their 
intuition, a more rational and methodical approach to dealing with this 
difficulty would be very useful. The probabilistic finite element method 
( PFEM) has been developed in response to these needs. 

In PFEM [1—3], finite element methods have been efficiently combined with 
second-order perturbation techniques to provide an effective method for 
informing the designer of the range of response which are likely in a given 
problem. The designer must provide as input the statistical character of the 
input variables, such as yield strength, load magnitude, and Young’s modulus, 
by specifying their mean values and their variances. The output then consists 
of the mean response and the variance in the response. Thus the designer is 
given a much broader picture of the predicted performance than with simply a 
single response curve. These methods are applicable to a wide class of 
problems, provided that the scale of randomness is not too large and the 
probabilistic density functions possess decaying tails. By incorporating the 
computational techniques we have developed in the past 3 years for efficiency, 
the probabilistic finite element methods are capable of handling large systems 
with many sources of uncertainties. 

Sample results for an elastic— plastic ten-bar structure and an elastic- 
plastic plane continuum with a circular hole subject to cyclic loadings with 
the yield stress on the random field are depicted in Figs. 1-4. For the ten- 
bar structure, a 5% coefficient of variation in the yield stresss gives a 13% 
coefficient of variation in the displacement of node 1 and an 11% coefficient 
of variation in the stress of element 1. For this example, along with many 
others (not shown here), PFEM compares very well with the Monte Carlo 
Simulation (MCS) and the Hermite Gauss Quadrature (HGQ) (see Fig. 2). This is 
an example where a situation where a small variance in the yield strength can 
result in a much larger variance in the response. It should be noted that the 
ratios of computer time are 1 to 400 when the PFEM is compared to MCS. 


*Work performed under NASA Grant NAG3- 535 administered by NASA Lewis 
Laboratories . 
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As for the elastic/plastic continuum problem, the mean displacement and 
stress are sinusoidal, resembling closely the forcing function (Fig, 4). The 
variances are close to zero until the plate begins to yield in compression. 
After this, the variance jumps to a higher value and remains steady until the 
yielding in tension begins. This phenomenon repeats every cycle. 

A third numerical example for PFEM methods is a turbine blade problem. 
Some results for the turbine blade model shown in Fig. 5 will be presented. 

The blade is subjected to a random impulsive load and the yield stress is 
random. 

A natural extension of these methods would be to consider fatigue and 
failure analysis. Finite element methods, such as PFEM, for analyzing fatigue 
and fracture in a probabilistic manner, are very scarce. The fracture related 
quantities such as fracture toughness, size and orientation of the cracks, are 
usually hard to determine exactly. These and other quantities, which govern 
the crack growth, can be treated by finite elements in a similar manner, 
although it woud be necessary to incorporate first and second order 
reliability methods and to embed singularities in the variational statements 
to correctly represent cracks. The experience obtained so far suggest that 
this is a logical extension of PFEM. 
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Problem Constants 

E = 3.0X1 0 7 lb/in 2 
Density = 0.3 lb/in 3 
Thickness = 1 .0 in 
L = 6.0 in 
R = 3.0 in 
Poissons Ratio = 0.3 
h y = 2000.0 lb/in 
w f = 1500.0 rad/sec 
Delta t = 1.0X10" 4 sec 
Rayleigh Damping Parameters 
e 0 = 0.0 e, = 1.5X10 -6 

Rand om Load 

24 Random Variables 
Coefficient of Variation = 0.1 
Mean Load = 2000.0 lb/in 
Spatial Correlation 
R(*i.Xj) = exp(-abs(x,— Xj)/Lp) 


M esh Data 

4 Node 2D Plane Strain Continuum 
Element in Radial Mesh 
784 Nodes, 720 Elements 
Point a = Element 1 
Point b = Element 15 
Point c = Element 346 
Point d = Element 360 
Point A = Node 1 
Point B = Node 16 
Point C = Node 385 
Point D = Node 400 

Rando m Ma t eria l 

15 Random Variables 
Coefficient of Variation = 0.1 
Mean Youngs Mod. = 3.0X1 0 7 lb/in 2 
Spatial Correlation 
R(*i,Xj) = exp(-abs(X|-Xj)/Le) 


Fig. 3 Problem Statement: Plain Strain Continuum with a 

Circular Hole. 


157 





MEAN STRESS (xIO 3 ) MEAN ( x10 ~ 3 ) 



Fig. 4a Mean and Variance of Node 400 y-Displacement versus 
Load Steps. 



Fig. 4b Mean and Variance of Stress in Element 15 versus 
Load Steps. 
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PROBABILISTIC SSME BLADES STRUCTURAL RESPONSE 
UNDER RANDCM PULSE LOADING 

Michael Shiao, Robert Rubinstein, and Vinod K. Nagpal 
Sverdrup Technology, Inc. 

Lewis Research Center Group 
Cleveland, Ohio 


The purpose of this work is to develop models of randcm inpacts on a Space 
Shuttle Main Engine (SSME) turbopunp blade and to predict the probabilistic 
structural response of the blade to these inpacts . The randan loading is 
caused by the inpact of debris. The probabilistic structural response is 
characterized by distribution functions for stress and displacements as 
functions of the loading parameters which determine the randan pulse model. 
These parameters include pulse arrival, anplitude, and location. The 
analysis can be extended to predict level crossing rates. This requires 
knowledge of the joint distribution of the response and its derivative. 

The model of randan iirpacts chosen allows the pulse arrivals, pulse 
arrplitudes, and pulse locations to be random. Specifically, the pulse 
arrivals are assumed to be governed by a Poisson process, which is 
characterized by a mean arrival rate. The pulse intensity is modelled as a 
normally distributed randcm variable with a zero mean chosen independently at 
each arrival. The standard deviation of the distribution is a measure of 
pulse intensity. Several different models were used for the pulse locations. 
For example, three points near the blade tip were chosen at which pulses were 
allowed to arrive with equal probability. Again, the locations were chosen 
independently at each arrival. 

The structural response was analyzed both by direct Monte Carlo simulation 
and by a semi -analytical method. In the Monte Carlo method, appropriate 
random number generators were used to develop simulated pulse arrival 
processes . These processes were used as forcing functions in a dynamic 
analysis of the SSME blade implemented by the computer code STAEBL 
(Structural Tailoring of Engine Blades). The dynamic analysis originally 
used by this program was a modal superposition based on up to five modes; 
this analysis was replaced by a direct time integration which used the 
Newmark Beta Algorithm. 

In the semi -analytical method, the classic analysis of shot noise by S. 0. 
Rice was generalized to the present problem. This analysis requires that the 
unit pulse response of the blade be known at each point where a pulse can 
arrive. The unit response was found by numerical simulation using the STAEBL 
code. Once the unit pulse response is known, all required distribution 
functions are developed analytically. 

Ccmparisons between the Monte Carlo studies and the semi-analytical method 
showed excellent agreement. Of course, the semi -analytical method has the 
advantage of requiring considerably less caiputer time to inplement. 
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An analysis was begun of the joint distribution of response and its time 
derivative. This joint distribution can be used to predict level crossing 
rates which are frequently used in fatigue life predictions. The joint 
distributions were evaluated both by direct simulation and sesni- 
analytically. In this case, the Monte Carlo simulations require very long 
integration times to produce smooth distributions. 

Possible extensions of this work include clustering effects in level 
crossings, direct simulation of the level crossing process, and analysis of 
the distribution of extreme values. 


OBJECTIVE 


RANDOM LOADINGS 


Estimate the influence of random loading on 
SSME blade responses 



1. Pulse 

2. Pressure 

3. Temperature 

4. Centrifugal 

RANDOM IMPULSE LOADING 

o Poisson arrival pattern 
o Random impulse amplitude 

o Random impulse location 

RANDOM STRUCTURAL RESPONSES 


SSME BLADE 


o Natural Frequency 
o Root Stress 
o Tip Displacement 
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RANDOM RESPONSES DUE TO RANDOM IMPULSE LOADING 


PROBABILITY CHARACTER OF SSME BLADE RESPONSE- 


SUBJECTED TO RANDOM PULSE LOADING 


• Three Response Models 


• Model 1 


Y(t) = i h(t-t k ) 


(1) Yj(t) 

h(t) 

t k 

(2) Y 2 (t) 

a k 


Z h(t-t k ) 

Unit Impulse Response 
A Poisson Process 

X a k h(t-t k ) 

Random Amplitude 


• Probability Density Function of Y(t) 
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Joint Probability Density Function of Y(t) and Y(t) 
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(3) Y 3 (t) 
h k (t) 


2 a k h k (t-t k ) 

Ranodm Unit Impulse Response 
(Random Pulse Location) 


Probabilistic characteristics of model 1 
and 2 are partially known 


where X 
T 

g k (y) 

g^(y.y) 


Poisson Arrival Rate 

Unit Impulse Response Duration 

k-Fold Convolution of Probability Density 

function of h(t) 

k-Fold Convolution of Joint Probability 
density function of h(t) and h(t) 
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QUANTIFYING UNCERTAINTIES IN THE STRUCTURAL RESPONSE OF SSME BLADES 


Vinod K. Nagpal 
Sverdrup Technology, Inc. 
Lewis Research Center Group 
Cleveland, Ohio 


A probabilistic study on turbopump blades of Space Shuttle Main Engine (SSME) 
undertaken at NASA Lewis is reaching the three year mark. The ultimate 
objective of this study is to evaluate the effects of random variations, 
generally called uncertainties, in the geometry, material properties, loading 
and their probabilistic combinations on the structural response of the blade. 
The results obtained on geometric and material properties uncertainties 
demonstrate that a methodology using probabilistic structural analysis 
methods appears to be a powerful and unique approach to quantify 
uncertainties . 

The application of this methodology can be extended to quantify the effects 
of random variations for other structural components also. Therefore, the 
development of the methodology has been discussed here in a rather 
generalized manner. This methodology consists of the following steps: 


( 1 ) The geometry and material properties can be randomly perturbed to 
simulate the realistic uncertainties. The simulation can be conducted using 
random numbers and using advanced perturbation techniques such as Monte 
Carlo's simulation. The randan numbers can be selected to have any given 
probabilistic distributions which have known statistical properties. In this 
study, the normal distribution was selected with known mean and standard 
deviation. 

(2) An analysis technique such as finite elements method can be used to 
estimate the structural response and provide means for discrete 
perturbations. The blade geometry and material properties were perturbed on 
the node and the element basis, respectively. 

(3) Statistical experiment designs such as full factional design can be used 
to determine the effects of study variables. The advantage of using this 
design is that it is cost effective and provides the estimates of the effects 
of individual study variables and their interactions. 

(4) The effects of all the study variables and their interactions can be 
evaluated for their significance using statistical tests. Further the 
probabilistic models can be developed to predict a mean response for given 
variations in the study variables. Probability distributions for response 
variables can also be developed for estimates of their range of variations. 
For this study, t-test, F-test, and^-^-test were used. 
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To quantify the uncertainties associated with the geometry and material 
properties of a SSME turbopump blade, a computer code known as STAEBL was 
used. A finite element model of the blade used 80 triangular shell elements 
with 55 nodes and five degrees of freedom per node. The whole study was 
simulated on the computer and no real experiments were conducted. The 
structural response has been evaluated in terms of three variables which are 
natural frequencies, root (maximum) stress, and blade tip displacements. 


The nodal coordinates (x, y, z) of the finite element mesh were perturbed to 
simulate the geometric uncertainties. The numbers of the material property 
matrix were perturbed for each element to simulate the material properties 
uncertainties. The perturbations were generated by a random number generator 
with preselected means and standard deviations. The magnitudes of means and 
Standard deviations for perturbations of both geometry and material 
properties, in this study, were taken as ten percent or less of the original 
values. These magnitudes were selected based on previous experimental 
results and experience. 

The results of the study indicate that only the geometric uncertainties have 
significant effects on the response. Uncertainties in material properties 
have insignificant effects. Also, the material properties interaction 
effects, which were created by variation in both material properties and 
geometry together have been found to have insignificant effects. A set of 
probabilistic models has been developed to predict the structural response 
for any given variations in geometry and material properties. Separate 
probabilistic models for only geometry variations have also been developed. 
Statistical tests indicate that these models are good fits. 
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OBJECTIVE 


CONDUCT PROBABILISTIC STRUCTURAL ANALYSIS 
OF SSME BLADES TO QUANTIFY UNCERTAINTIES 
ASSOCIATED WITH: 

GEOMETRY 

MATERIAL PROPERTIES 


WHY? 

EVALUATE TOLERANCE LIMITS 

- HIGHER TOLERANCE 
LOWER COST 


BLADE GEOMETRY 

SSME HIGH-PRESSURE FUEL TURBOPUMP 
1st STAGE TURBINE BLADE 



METHODOLOGY 

• Selected random distribution to simulate 
uncertainties 

Normal 

• Modelled blade with finite element model 

80 Elements 
55 Nodes 

• Used experiment design to perform simulation 

Full factorial design 

• Analyzed response 

Probabll Istic Models 
Probabilistic Distributions 
Statistical Tests 

Figure 2. 


Figure 1. 
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RANDOM DISTRIBUTIONS 



FINITE ELEMENT MODEL OF SSME BLADE 
Figure 3 . 
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OPTIMIZED HAY OF STUDYING THREE VARIABLE EFFECTS 
FULL FACTORIAL DESIGN 

Figure 4. 


STRUCTURAL RESPONSE 

• NATURAL FREQUENCIES 

First 

Second 

.Third 

• ROOT (MAXIMUM) STRESS 

• TIP DISPLACEMENTS 


Figure 5. 
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TABLE 1 - PROBABILISTIC MODELS - GEOMETRIC P ERTURBA T I ONS 
MODEL: 

OEP. VAR. • CONSTANT ♦ CO OT ♦ COOT ♦ COOT * COOT » * COOT *0^ ♦ COOT 


DEPENDENT 

VARIABLE 

CONSTANT 

COEFFICIENTS OF 





<*4 

<?5 

^6 

FIRST FREQ. 

6105.0 

-032 

-2915 

-607 

-8362 

1897 

-9348 

SECONO FREQ. 

9475.1 

-1374 

-2961 

-400 

-11398 

-3869 

-9095 

THIRD FREQ. 

15792. 

-9434 

27592 

-7944 

-17663 

-4899 

-44147 

ROOT STRESS 

63323. 

49707 

103320 

4177 

68960 

-48191 

2B3960 

TIP OISPL. 

.00196 

.0119 

. 1267 

.0206 

-.0132 

-.0419 

.0067 


F tests indicated that all models are good fits. 


DISTRIBUTION OF NATURAL FREQUENCIES 
- MATERIAL PROPERTIES PERTURBATIONS 



FIRST NATURAL FREQUENCY 


Figure 6. 


172 









STATISTICAL TESTS 


• t-Test 

t = Value - Mean 
St. Dev. 

• F-Test 

F = (Residual Sun of Squares)/p.F . 
(Regression Sum of $quares)/D.F7 

I Qc 2 Test 

'V 2 = (Value - Mean) 2 
** — Variance — 

l Plots 


Figure 7 . 


CONCLUSIONS 


• Methodology for SSME Blade has been developed and 
appl led. 

• Methodology can be extended to other structural 
components. 

I Geometric uncertainties showed significant effect. 

• Material properties uncertainties and their Interactions 
have Insignificant effects. 

• Range of variation In response quantified. 

Figure 8. 
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COMPOSITE LOAD SPECTRA 

FOR SELECT SPACE PROPULSION STRUCTURAL COMPONENTS* 

J . F . Newe 1 1 

Rocketdyne Division, Rockwell International 
Canoga Park, California 


The objective of this program is to develop generic load models to simulate 
the composite load spectra (CLS) that are induced in space propulsion system 
components representative of the space shuttle main engines (SSME). These 
models are being developed through describing individual component loads with 
an appropriate mix of deterministic and state-of-the-art probabilistic models 
that are related to key generic variables. Combinations of the individual 
loads are used to synthesize the composite loads spectra. 

A second approach for developing the composite loads spectra load model 
simulation, the option portion of the contract, will develop coupled models 
which combine the individual load models. Statistically varying coefficients 
of the physical models will be used to obtain the composite load spectra. 

The need for this type of technology advancement is apparent from the demand 
for higher performance, lighter weight, components that yield higher operating 
pressures, temperatures, vibration and flow loads. Figure 1. The difficulty 
in installation, cost, and potential for new failure mechanisms limit the 

required instrumentation to adequately define or verify loads. Proper 
quantification of loads are required to minimize adding conservatism on 
conservatism and to specify results so they can be accurately used for 

structural analysis and risk assessment. 

Rocketdyne is teamed with Battelle to develop this methodology. Rocketdyne is 
using their background and expertise for the advancement in the loads 
definition and analysis methods along with the extensive SSME test database 

for comparitive information. Battelle is utilizing their expertise to develop 
an advanced probabilistic code to represent the individual loads and load 
spectra. Rocketdyne is packaging this total work in an easy to use expert 
system type code. The following discussion covers an overview and 

Rocketdyne's effort. Battelle's contribution is presented in a separate 

report. 

The CLS development is a 3-year base program with a 2-year option program. 
The first 2 years of the base program have been completed. The effort has 
three major tasks: probabilistic model theory and development, code 

development and code validation and verification. Figure 2. Four classes of 
rocket engine components are being used as examples for the load 
development — turbine blades, transfer ducts, LOX posts and engine system 
ducts. Figure 3 and 4. The available SSME instrumentation for use in 

developing SSME related statistical data and load verification are also shown 
in Figure 3. 


*work performed under NASA Contract NAS3-24382 

1 7 5 


PRECEDING PAGE BLANK NOT FILMED 


The probabilistic modeling requires simulation of the individual loads 
summarized in Figure 5. Several shape simulations are required to 
characterize these loads. The start and cutoff transients include a nominal, 
random over nominal and spike load shape. The overall duty cycle required 
to define engine operating conditions is represented as a nominal shape with a 
random variation about the mean value. The basic nominal shapes of the duty 
cycle operational requirements and inlet parameter operation bounds are 
controlled by the engine contract requirements. These requirements are based 
on vehicle thrust and operation needs. Figure 6. Environments like vibration 
have a nominal variation dependent on a components power and/or speed, random 
variation from hardware geometry and test to test conditions and transient 
spikes from side loads, pops and chugs. The steady state dynamic loads are 
represented in the frequency domain for both random and sinusoidal loads. 
Figure 7. 

Since the CLS methodology is to be applied to advanced engines, the loads are 
being developed from a generic basis. The use of hardware specific random 
variations like pump efficiency, head rise or flow resistance allow for engine 
to engine load variations. Accounting for variations in power level, pump 
speed, and engine inlet operating conditions allow for test to test variation. 
Figure 8. The use of key variables and appropriate physical and probabilistic 
models assures the methodology is applicable to variations in the SSME or to 
new rocket engine designs. The results obtained from the loads spectra models 
are being compared with available SSME engine measured results and available 
analytical calculations. 

The probabilistic loads model being developed by Battelle is implemented as 
part of an expert system developed as part of the program. The expert system 
is a tool to generate and analyze composite loads of a rocket engine design 
and to supply these loads for use in either deterministic or probabilistic 
finite element computer codes to perform structural analysis of engine 
components. The probabilistic models are generic; the statistical information 
utilized is primarily from the SSME test database. Expert opinion and other 
engine background data are used where appropriate to complete the loads 
picture. 

The knowledge-based system manages the database, provides expert knowledge 
relative to the generic probability loadings and generates the individual and 
composite loads, Figure 9. The Battelle developed ANLOAD module performs the 
probabilistic modeling and statistical analysis. A database system has been 
developed to efficiently represent the knowledge. This database system 
facilitates the communication between the expert system and the knowledge 
base. The current knowledge base contains information for SSME type engines, 
system type loads using the influence coefficient method, and local turbine 
blade load scaling methods. 

The expert system is a rule based system. The rules are modularized where 
each module is designed to solve a particular problem or to perform a task. 
The load expert system LDEXPT Version 2.0, Figure 10, has rule modules for 
engine system dependent loads for all four components as well as selected 
individual local component loads. The rules so far mostly relate to overall 
process control and information retrieval. Figure 11. The on-going rule 
development is working on local individual load components and the more 
complex composite load spectra. Knowledge for the transfer ducts is partially 
developed and is being added to the system. The other two component loads are 
being developed. 
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FIGURE 2 

COMPOSITE LOADS SPECTRA 
BASE PROGRAM 
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COMPONENT IMPLEMENTATION AND 

LOAD COMPLEXITY IMPROVED IN EACH CODE VERSION 


FIGURE 3 

CLS ANALYSIS 

COMPONENTS AND STANDARD INSTRUMENTATION 


ORIGINAL PAGE IS 
OF POOR QUALITY 
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FIGURE 4 

HIGH PRESSURE PUMP DISCHARGE DUCT 


ORIGINAL PAGE IS 
Of* POOR QUALITY 
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SUMMARY MATRIX OF INDIVIDUAL LOADS VS COMPONENTS 
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LOW FREQ. & TRANSIENT 





MECHANICAL VIBRATION 
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FREQUENCIES DEPENDENT ON PUMP SPEED 





FIGURE 9 

LDEXPT LOAD EXPERT SYSTEM 
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IMPLEMENTING PROBLEM-SOLVING KNOWLEDGE 
• E.G. SELECTING INDEPENDENT LOADS BASED ON GAINS 


LDEXPT-LOAD EXPERT SYSTEM 



86 


PROBLEM TEXT FILE LOAD DATABASE 

RULE FILE 














FIGURE 11 

LDEXPT LOAD EXPERT SYSTEM 
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PROBABILISTIC LOAD MODEL DEVELOPMENT AND VALIDATION FOR COMPOSITE 
LOAD SPECTRA FOR SELECT SPACE PROPULSION ENGINES 

R. Kurth 

Battelle Memorial Institute 
Columbus, Ohio 

and 

J.F. Newell 

Rocketdyne Division, Rockwell International 
Canoga Park, California 


A major task of the program to develop an expert system to predict the 
loads on selected components of a generic space propulsion engine is the 
design, development, and application of a probabilistic loads model. 

This model is being developed in order to account for the random nature 
of the loads and assess the variable load ranges’ effect on the engine 
performance. 

There are several requirements of the probabilistic load model which 
restrict the use of some probabilistic methodologies. First, the 
methodology must be capable of addressing four different shapes for the 
sample paths in the stochastic processes describing the loads: nominal, 
random, spike (rare event), and periodic. This implies that the model 
must be capable of approximating a variety of different forms for the 
load processes, some, or all, of which may be non- stationary. Second, 
the model must be applicable to generic engines operating under a wide 
variety of possible mission history profiles, typified by the demanded 
power levels. This requirement implies that the model must be capable of 
dealing with both random variables, as well as stochastic processes, that 
may not be describable by standard distributional forms, for example 
bimodal behavior or discrete (on-off) processes. Also, there are sixty 
to seventy individual loads and/or engine parameters that can affect the 
composite load calculations (although, usually only ten to thirty of 
these are important for a single, specific composite load calculation), 
and, therefore, the methodology must be as efficient as possible without 
incurring a substantial loss of accuracy. Finally, the methodology is 
required to provide varying levels of accuracy, the level to be user 
specified. 

A probabilistic model capable of addressing each of these topics has been 
developed. The model is based primarily on simulation methods, but also 
has a Gaussian algebra method (if all variables are near normal), a fast 
probability integrator routine (for the calculation of low probability 
events), and a separate, stand alone program for performing barrier 
crossing calculations. Each of these probabilistic methods has been 
verified with theoretical calculations using assumed distributional 
forms. The results of the these verification studies are discussed in 
the full presentation. 
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The governing philosophy for the probabilistic methodology is that each 
process can be described by a probability density function. To implement 
this philosophy for stochastic processes, approximate methods for 
modeling the mission history profile are needed during portions of the 
mission. Independent (controller), individual loads and engine 
parameters are treated as random variables whose distributional form and 
parameters are determined by the controller demanded power level 
(thrust). The dependent, individual loads and composite loads are then 
modeled by assuming a quasi-steady state system. This implies that 
during the power up and down, and throttle up and down phases of a 
mission the time step size be sufficiently small to warrant a steady 
state approximation, from the probabilistic point of vi ew. That is, the 
physics may not be in a quasi-steady state, but it can be assumed that 
the loads of interest are stationary during this time. Between t ime 
steps the non- stationary loads are adjusted according to the data base 
information, or user supplied input. An example of the coupling of 
quasi-steady state and steady state phases will be discussed. 

For transient events, such as engine start and cut-off, the changes in 
the magnitude of the loads is both large (relatively) and of short 
duration (one to four seconds). For these events the load profile is not 
assumed to be quasi-steady in order to keep the time step size relatively 
large and reduce the computational time. Instead, the fluctuations about 
the nominal response are ignored and the random nature of the transient 
is contained in the variable amplitude and time of occurrence. The load 
is modeled as a piecewise linear response over the duration of the 
transient. For example, in a two-stage engine, the temperature in the 
fuel pump is expected to have two or three spike values before settling 
down to a relatively steady rate of increase to the full power operation 
temperature level. The occurrence, or non-occurrence, of the third peak 
is handled statistically as a random variable, as well as the amplitude 
of the peak and the time at which the peak amplitude occurs. Examples of 
the transient model, and its coupling to other mission phase types is 
provided. 

Because the probabilistic model is treating the problem in the time 
domain, special attention must be paid to periodic, or sinusoidal, loads. 
Since their magnitudes are usually more sensitive to frequency than they 
are to any time dependency (other than the time dependency of the 
frequencies), they are treated separately in the probabilistic model. 

At this time most of the available data for the model validation has been 
obtained from the space shuttle main engine program. The data which has 
been examined thus far appears to agree, within statistically accuracy, 
with the model predictions obtained from the validation studies. Further 
testing with other engine data is planned to demonstrate the generic 
capability of the model. An example of a validation analysis is 
provided . 
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INTRODUCTION 


Probabilistic Load Analysis 

- The primary goal is the development of a probabilistic 
methodology that is capable of handling any generic space 
propulsion engine 

— Requires that any anticipated distribution of engine load or 
parameters must be addressable by the model 

- Must be able to analyze, in a manner consistent with the 
stochastic process being modeled, either linear or non-l(near 
ombinations of the inputs which may be non_ stationary 

Probabilistic Load Model 

The above points have led to the development of a model based 

on the following points: 

• Tabulation of means, variances, and distribution types for 

• Scaling laws, or expert opinion, for modifications of data 
base for engines currently not included 

• A probabilistic model that provides varying levels of 
sophistication for calculating the stochastic load processes 

• Capability to analyze nominal, random, or periodic loads 


PROBABILISTIC LOAD MODEL 


Stochastic Processes 

— The stochastic load processes are modeled using PDPs to 
describe the individual load variables and engine parameters 

— There are four primary shapes that are modeled: 

(a) Nominal 

(b) Random 

(c) Spike 

(d) Peridoic 

- Because the load is a stochastic process there are three 
different methods for displaying the results: 

(a) Nominal plot 

(b) Sample function plot 

(c) Probability tables of exceeding specified load levels 
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PROBABILISTIC LOAD MODEL 

(continued) 


Individual loads 


- Probability density functions 

— Data analysis, or calculated using engine influence coefficient 


Composite loads 


— User supplied: 

Mission history profile 
Accuracy required 
Output desired 

- Composite load is calculated from the probabilistic method 
selected from the input accuracy needed and the desired output 


Non-stationary processes 

— Computer program decides if the individual loads are non- 
stationary (This option can be suppressed). If they are not 
stationary, then new distribution parameters are calculated 
based on the current power level 

- Relatively small time step is used to approximate a stationary 
process over the time increment 


Probabilistic Methodology 


— Methodology based primarily on simulation methods 

— Gaussian algebra, fast probability integrator (Wirsching and 
Wu), and barrier crossing methods are also available for 
specialized calculations 

— Distribution fitting subroutine is included in the model to 
provide summary statistics on the stochastic processes or 
random variables 


Displayed Output 


- Discrete load level, with its associated probability 

- Updated distributions for the individual load (where warranted) 
and composite loads at specified times 

- Interface with printer plot or LOTUS plotting packages 

- All output is stored on a file called DPD.OUT 


Data Base Update 


— If data analysis for individual load has been performed, it is 
inserted in the data base (not currently automated) 

— The output for the composite load will be cataloged so that 
future users may simply refer to this analysis without having 
to re-do the entire calculation 
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Probabilistic Model Flowchart 
Model Input and Set-up 
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Probabilistic Model Flowchart 
Modei Calculations 
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Sample Calculation i 

RASCAL Verification 


Prediction of Sum Of 3 Normal Variables 


RASCAL Prediction 



— Sum of three normal variables N(10 f 1) 

— Known mean of 30 and variance of 3 

— Prediction obtained from the 1% to 99% values of the CDF is 
within 0.02% relative error of the theoretical distribution 
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Cumulative Diatribution Function 


Sample Calculation I! 

Low Probability Calculations 

RASCAL Prediction: Failure Probability 


Exceedance Of Yield Strength 



Yield Strength - Von Mises Stress 

+ 10-100 o 40-25 x 250-40 Chin -Lind 


— Pressure vessel failure probability model used because previous 
analysis is available for comparison 

— Fast probability integerator (Wirsching and Wu) 

— RASCAL and FPI predictions closely agree 
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HCPFTP Temperature 
(Thousands ) 


ORIGINAL PAGE i3 
OF POOR QUALITY 


SGimple Calculation II! 

Transient Load Model 

Transient Temeperature Model 



□ Transient Model + Average Bounds 


- Frequency of: 

Number of peak amplitudes 
Magnitude of peak amplitudes 
Timing of peak amplitudes 

- Piecewise linear model is used 

- Plot compares model to ensemble average for the verification. 
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H PFTP Temperature 
(Thousands) 


Sample Calculation IV 

Transient Model Validation 


Transient Temeperature Model 



Tlmt After Engine (seconds) 

□ 902349 f 902356 & 902363 


— SSME HPFTP data was examined for data sets that spanned the 
the possible range of sample paths 

— Plot compares transient model and three selected test sets 

— Comparison appears very good 
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Sample Calculation V 

Mission Phase Linking 


ORIGINAL 

of poor 


PAGE 13 

QUALITY 


Transient Temeperature Model 



— Transient spike values not shown 

— Influence function only valid above 65% rated power level 

— Extrapolation used to move from the end of the tranmsient 
data set to the beginning of the quasi— steady phase 

— Agreement between the point reached from extrapolation and 
the 65% level for temperature is very good 


SUMMARY 

• A probabilistic load model has been developed that is 
capable of addressing all portions of the mission profile 

• The model can address constant, random, spike, and 
periodic load forms 

• Comparisons to date have shown the model to work well 
in comparison to available data 
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STRUCTURAL TAILORING USING THE SSME/STAEBL CODE* 


Robert Rubinstein 
Sverdrup Technology, Inc. 
Lewis Research Center Group 
Cleveland, Ohio 


Space Shuttle Main Engine (SSME) blades are subject to severe thermal, 
pressure, and forced vibration environments. An SSME blade design must meet 
tight clearance, fatigue life, and stress limit constraints. Because of the 
large number of potentially conflicting constraints, a "manual" design 
P ro F e 5* ur ® ^ require many time-consuming iterations. Structural 
optimization provides an automated alternative. Any number of analyses, 
design variables, and constraints can be incorporated in a structural 
optimization computer code. This idea has been applied to develop the code 
SSME/STAEBL, which is a stand-alone code suitable for automated design of 
SSME turbopump blades. 

SSME/STAEBL was developed by systematically modifying and enhancing the 
STAEBL (Structural Tailoring of Engine Blades) code developed by Pratt and 
Whitney under contract to NASA Lewis Research Center. STAEBL was designed 
for application to gas turbine blade design. Typical design variables 
include blade thickness distribution and root chord. Typical constraints 
include resonance margins, root stress, and root to chord ratios. In this 
program, the blade is loaded by centrifugal forces only. 

Additions and modifications of STAEBL included in SSME/STAEBL include (1) 
thermal stress analysis, (2) gas dynamic (pressure) loads, (3) temperature 
dependent material and thermal properties, (4) forced vibrations, (5) tip 
displacement constraints, (6) single crystal material analysis, (7) blade 
cross section stacking offsets, (8) direct time integration algorithm for 
transient dynamic response. Capabilities are also included which permit data 
transfer from finite element models and stand-alone analysis. 

Several design optimization studies have been completed using an SSME blade 
design to test these various capabilities. Optimization studies have been 
completed to test the influence of thermal and pressure loads and temperature 
dependent properties on optimal blade design. Comparison between designs 
optimized under centrifugal loads only and under centrifugal, thermal, and 
pressure loads with temperature dependent blade properties shows that the 
additional loads require additional weight to meet all design constraints. 
The difference between the designs can be attributed to material property 
temperature dependence, which in this case forces a much tighter root stress 
constraint . 

Design optimization studies for a blade made of a typical single crystal 
material showed relatively little effect of crystal axis orientation. This 


*Work performed under NASA Lewis Task Order No. 5207. 
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study was dominated by a root stress constraint which was violated by the 
initial design. It was found that root stress is influenced much less by 
crystal orientation than by the geometric design variables. The result is 
that the optimized design is found by adjusting the blade geometry 
significantly, but the crystal axis orientation insignificantly. Of course, 
in blade designs dominated by natural frequency constraints in particular, a 
different conclusion could obtain. 

Optimization studies were undertaken to assess the influence of blade cross- 
section offsets on optimal design. These offsets are defined as the 
differences between the blade cross section centers of gravity and a straight 
line perpendicular to the engine axis. Blade designers use these variables 
to balance centrifugal and pressure loads. In a typical optimization study, 
the initial design violated the root stress constraint. Root stress is 
extremely sensitive to offset variables, and the optimized design recommended 
a significant change in stacking. 

It can be concluded that structural optimization provides a practical and 
flexible approach to SSME blade design. The general structure of 
optimization algorithms allows many design variables, constraints, loads, and 
analysis options to be incorporated in the optimization procedure. In the 
design optimization studies reported here, which were dominated mainly by 
root stress constraints, it was found that temperature dependent properties 
can require additional weight in the optimized design, mainly by reducing the 
yield stress. Crystal axis orientation was relatively unimportant design 
variable in these studies. Blade cross-section stacking offsets were shown 
to have very significant effect on root stress-dominated designs. This 
conclusion agrees with the general experience of blade designers. 

OBJECTIVE 


Apply structural optimization methods to 
SSME blade design. 


smuciim/u. lAiionmo or engine blades isiaebd 
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OPTIMIZATION STUDY - SSME TURBOPUMP BLADE 
COMPARISON BETWEEN INITIAL AND OPTIMIZED DESIGN 


Initial design 


Optimized design under 
centrifugal loads only 

Xspan THK (in) 

chd (in) 

X span 

THK (IN) CHD (IN) 

0. .233 
50. .138 
100. ,092 

.761 

i: 

.22<i .890 
.082 .681 
,065 .650 

Natural frequencies (cps) 

3562. 

R583. 

7916. 


3954. 

48G8. 

7969. 

Root stress (ksi) 

S3. 


108. 

Tip displacements 

Untwist (deg) 
Uncamber (deg) 
tip Ext (in) 

iil! 


0.6&27 

Blade Weight (lb) 
Number of blades % 
Stage weight (lb) 

j 




OPTIMIZATION STUDY - SSME TURBOPUMP BLADE 
EFFECT OF THERMAL, PRESSURE LOADS AND FORCED VIBRATION ON OPTIMAL DESIGN 


Centrifugal loads only, Representative thermal, pressure loaos 

TEMPERATURE INDEPENDENT PROPERTIES TEMPERATURE DEPENDENT PROPERTIES 


X SPAN 

0 . 

50. 

100 . 


THK (IN) 

.22R 

.082 

.065 


CHD (IN) 

.890 

in 


Natural Frequencies (cm) 

88 : 

7969. 

Root Stress (ksi) 108. 


X span 

THK (IN) 

chd ( 

0. 

.228 

.89 

50. 

.082 

.68 

100. 

.077 

.65 


8835. 

65. 


Tip Displacements 


Untwist (deg) 
Uiicamber (deg) 
Tip Ext (ini 



0.0i09 


Forced Response 
Margins 



:82 

.010 


Blade weight (is) 
Number of Blade s. 
Stage weight (lb) 


.093 


J.18 
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SINGLE CRYSIAL ANALYSIS 


EFFECT OF CRYSTAL ORIENTATION ON OPTIMAL DESIGN 


Crystal aligned Crystal aligned 95* 


WITH 

ROOT CHORD 


TO ROOT 

CHORD 

X SPAN 

THK . (IN) 

CHD . (IN) 

THK. (IN) 

CHD. (IN) 

0 

.23 

.89 

.23 

.89 

50 

.11 

.69 

.08 

.69 

100 

.06 

.65 

.06 

.65 


frequencies 

(cps) 5110. 

6851. 

13590, 


5079. 

6897. 

13870. 


ROOT STRESS 

( KS 1 ) 63. 


65. 


STAGE HEIGHT 
(LB) 


3.26 


3.15 


Crystal orientation Euler angles 

d .89 E-3 95.00 

(5 .57 E-2 -.20 E-l 

Y -.51 E-3 .98 E-2 


SINGLE CRYSTAL ANALYSIS 

Gradient of root stress constraint at first iteration 


VARIABLE 


GRADIENT 

THICKNESS 

1 

-.29 

THICKNESS 

2 

-.16 

THICKNESS 

3 

.01 

THICKNESS 

9 

.09 

THICKNESS 

5 

.01 

ROOT CHORD 

1.55 

EULER ANGLES 




.02 

P 


-.01 

Y 


.00 
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Optimize! Ion Study Including blode stocking dfslgn vorlobles 



lull lal 

Design 

Final 

Design 


X spoil 

Ilik 

did 

X soon 

tl«k 

ClKf 

0 . 

.233 

l.O'iI 

0 . 

.276 

.890 

'50. 

.1)8 

.80*1 

50. 

.1)3 

.688 

100. 

.092 

.76) 

100. 

.091 

.650 

blode stocking 





A - 0. 

B ■ 

0 . 

A « 

-.0012 

i - .0078 

C - 0. 

0 • 

0 . 

C - 

.00)5 

D - .0095 


Root Stress 

82* ksl 7* kif. 


2 ox Is along soon# blade lengtli • l 

Ax • x displacement fra* center of gravity 

Ay * y displacement (ran center of grovlty 

Ax* A (l/U ♦ B It/It 1 
A i - C (I/O ♦ 0 (*/l» 2 


Comparison between eflecti of thickness. rool cliord. ond blode 
docking on root (lied. 


Optimization run with rool stress constraint vloloied by Initial detlgn 


Design Vorloble Crodleni of Bool Sired 


I hk 1 IrooO -•?! 

Ilik 2 

I lik J .00 

thk M *02 

Ilik S lllpl .00 


root chord 

A 9.20 

I B.22 

blode docking 2.39 

D 2.B2 
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NONISOTHERMAL ELASTO-VISCO-PLASTIC RESPONSE OF SHELL-TYPE STRUCTURES* 


G.J. Simitses, R.L. Carlson, and R. Riff 
Georgia Institute of Technology 
Atlanta, Georgia 


The prediction of inelastic behavior of metallic materials at 
elevated temperatures has increased in importance in recent years 
The operating conditions within the hot section of a rocket motor 
or a modern gas turbine engine present an extremely harsh thermo- 
mechanical environment. Large thermal transients are induced each 
time the engine is started or shut down. Additional thermal tran- 
sient from an elevated ambient, occur whenever the engine power 
level is adjusted to meet flight requirements. The structural 
elements employed to construct such hot sections, as well as any 
engine components located therein, must be capable of withstand- 
ing such extreme conditions. Failure of a component would, due to 
the critical nature of the hot section, lead to an immediate and 
catastrophic loss in power and thus cannot be tolerated. Conse- 
quently, assuring satisfactory long term performance for such 
components is a major concern for the designer. 

Traditionally, this requirement for long term durability has 
been a more significant concern for gas turbine engines than for 
rocket motors. However, with the advent of reusable space vehi- 
cles, such as the Space Shuttle, the requirement to accurately 
predict future performance, following repeated elevated tempera- 
ture operation, must now be extended to include the more extreme 
rocket motor application. 

A mathematical model and solution methodologies for analyzing 
structural response of thin, metallic shell-type structures under 
large transient, cyclic or static thermomechanical loads have 
been developed. Among the system responses, which are associated 
with these load conditions, are thermal buckling and creep buck- 
ling. Thus, geometric as well as material-type nonlinearities (of 
high order) can be anticipated and have been considered in the 
development of the mathematical model. Furthermore, this was 
accommodated in the solution procedures. 

A complete true ab-inito rate theory of kinematics and kine- 
tics for continuum and curved thin structures, without any res- 
triction on the magnitude of the strains or the deformation, was 
formulated. The time dependence and large strain behavior are 
incorporated through the introduction of th^ time rates of the 
metric and curvature in two coordinate systems; a fixed (spatial) 
one and a convected (material) coordinate system. The relations 
between the time derivative and the covariant derivatives (gra- 


*This work was performed under NASA Grant No. NAG3-53^i. 
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dients) have been developed for curved space and motion, so that 
the velocity components supply the connection between the equa- 
tions of motion and the time rate of change of the metric and 
curvature tenors. 

The metric tensor (time rate of change) in the convected ma- 
terial coordinate system is linearly decomposed into elastic and 
plastic parts. In this formulation, a yield function is assumed, 
which is dependent on the rate of change of stress, metric, tem- 
perature, and a set of internal variables. Moreover, a hypoelas- 
tic law was chosen to describe the thermoelastic part of the de- 
formation. 

A time and temperature dependent viscoplastic model was for- 
mulated in this convected material system to account for finite 
strains and rotations. The history and temperature dependence 
were incorporated through the introduction of internal variables. 
The choice of these variables, as well as their evolution, was 
motivated by phenomenological thermodynamic considerations. The 
nonisothermal . elastic-viscoplastic deformation process was des- 
cribed completely by "thermodynamic state" equations. 

The procedure employed permits the rates of the field for- 
mulation to be. interpreted as increments in the numerical solu- 
tion. This is particulary convenient for the construction of a 
finite element models together with incremental boundary condi- 
tion histories. Finite element formulation was developed for cur- 
ved beams and shells. The element matrixes were derived directly 
from the incrementally formulated equations using tensor oriented 
procedure. 

Finite element solution of any boundary-value problem invol- 
ves the solution of the equilibrium equation (global) together 
with the constitutive equation (local) . Both equations are solved 
simultaneously in a step by step manner. The incremental form of 
the global and the local equation can be achieved by taking the 
integration over the incremental time step bt=tj +1 -tj. The rec- 
tangular rule has been applied to execute the resulting time in- 
tegration. 

For structures with unstable deformation paths (snap-buckling 
phenomenon) , accurate and efficient description of the motion of 
the structure was obtained by inclusion of the inertia forces. 

Applications: The response of a clamped circular arch and of 
a cylindrical panel were studied. The shallow circular clamped 
arch subjected to a single central concentrated load, as shown in 
Fig. 1, was analyzed. The material chosen for the numerical expe- 
rimentation is the carbon steel C-45 (DIN 1720) with E = 10 7 psi, 
V= 0.3 and <i y = 2.7 *10 4 psi at room temperature. 

The arcn response, the deflection time history and the influ- 
ence of temperature on the arch response are shown in Figs. 2,3 
and 4, respectively. 

A thin cylindrical shell panel simply supported on all sides, 
made of the same material as the arch, and subjected to in-plane 
loads along the generators as shown in Fig. 5 was also studied. 

The deflection time history of the panel is shown in Fig 6, 
for a value of N|»|* = 20 lbs/in. This load is well below the criti- 
cal (buckling) load for the geometry which is 42.15 lbs/in. 
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FIG. 4 THE INFLUENCE OF TEMPERATURE 




5 CYLINDRICAL PANEL 


FIG. 6 CYLINDER PANEL CREEP RESPONSE 
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DYNAMIC CHARACTERISTICS OF SINGLE CRYSTAL SSME BLADES 

L.A. Moss and T.E. Smith 
Sverdrup Technology, Inc. 

Middleburg Heights, Ohio 


The SSME High Pressure Fuel Turbopump (HPFTP) blades are 
currently manufactured using a Directionally Solidified (DS) 
material, MAR-M-246+Hf . However, a necessity to reduce the 
occurrence of fatigue cracking within the DS blades has lead to 
an interest in the use of a Single Crystal (SC) material, PWA- 
1480. 

A study was initiated to determine the dynamic characteristics of 
the HPFTP blades made of SC material and find possible critical 
engine order excitations. This study examined both the first and 
second stage drive turbine blades of the HPFTP. The dynamic 
characterization was done analytically as well as experimentally. 
The analytical study examined the SC first stage HPFTP blade 
dynamic characteristics under typical operating conditions. The 
blades were analyzed using MSC/NASTRAN and the finite element 
model shown in Figure 1. Two operating conditions, 27500 RPM and 
35000 RPM, were investigated. Each speed had a unique blade 
temperature distribution. Bench experiments were conducted to 
determine the non-rotating natural frequencies and mode shapes. 

Parametric studies were done varying the crystal orientation with 
respect to the blade. Crystal orientation was important because 
of the anisotropic nature of the material. Youngs Modulus was 
greatest along the crystal’s < 1 1 1 > direction. The crystal <111> 
direction was referenced within the blade according to its 
projected rotate angle and tilt angle shown in Figure 2. The SC 
blades were initially ana. lyzed at a tilt angle of 54.74 degrees 
such that its crystal growth direction was aligned in the blades’ 
span direction, similar to the DS blade. 

For the first three modes the frequency appeared to be least 
dependent on the crystal rotation angle. For the same modes the 
crystal tilt angle had a greater influence on the modal 
frequencies as shown in Figure 3. For modes four through six the 
crystal rotate angle had a more significant effect on natural 
frequencies as shown in Figure 4. 

Campbell diagrams were constructed with the following critical 
engine orders determined by Rocketdyne; 11, 13, 15, 24, 26, 28, 
and 41. The modal analysis results for the SC blades were found 
to be worse than the DS when considering the engine order 
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interferences (see Figures 5 and 6 ). These results indicated 
critical engine order crossings may occur over the first three 
modes of the SC blades. 

Additional investigations attempted to determine an optimum 
crystal orientation which would most effectively avoid the above 
interferences. The procedure was to vary the orientation of the 
crystals <111 > direction between the extremes. The extremes 
would be strongest crystal direction either blade spanwise or 
chordwise. The best orientation when considering the a +/- 5% 
deviation of excitation frequency occurred at a. tilt angle of 90 
degrees and a rotate angle of 0 degrees (see Figure / ). It 
appears that for any orientation, the third mode interference 
will exist ( Figure 8 ). 


MSC/N ASTR AN FINITE ELEMENT MODEL 



MATERIAL ORIENTATION OF SINGLE CRYSTAL BLADE 



CRYSTAL. GROWTH 
DIRECTION <001> 



FIGURE 2 
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SSME HPFTP - 1st Stage Blade 
PWA-1400; 35000 RPH: Taf-1400 F. and Tsh-1200 F. 
1st Mode 


Frequancy (Hz) 



FIGURE 3 


SSME HPFTP - 1st Stage Blade 
PWA-1480; Nonrotating; 70 degress F. 
4th Mode 

Frequency (Hz) 



FIGURE 4 


SSME HPFTP - 1ST STAGE BLADE 
CAMPBELL DIAGRAM; DS MAR-M-246 Hf MATERIAL 



FIGURE 5 
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SSME HPFTP - 1ST STAGE BLADE 
CAMPBELL DIAGRAM: PWA-1480 MATERIAL 
SC ORIENTATION: TILT-54.74. ROTATE-O.O 

FREQUENCY (Hz) 



FIGURE 6 


SSME HPFTP - 1ST STAGE BLADE 
CAMPBELL DIAGRAM: PWA-1480 MATERIAL 
SC ORIENTATION: TILT«90.00. ROTATE-OO.O 

FREQUENCY (Hz) 



FIGURE 7 


SSME HPFTP - 1st Stage Blade 
PWA-1480: 35000 RPM: Taf-1400 F. and Tsh-1200 F. 
3rd Mode 



FIGURE 8 
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SSME BLADE DAMPER TECHNOLOGY 

Robert E. Klelb 
NASA Lewis Research Center 
Cleveland, Ohio 

and 

Jerry H. Griffin 
Carnegie-Mellon University 
Pittsburgh, Pennsylvania 


Before 1975 turbine blade damper designs were based on experience and very 
simple mathematical models. Failure of the dampers to perform as expected 
showed the need to gain a better understanding of the physical mechanism of 
friction dampers. Over the last 10 years research on friction dampers for 
aeronautical propulsion systems has resulted In methods to optimize damper 
designs . 

The first-stage turbine blades on the SSME high-pressure oxygen pump have 
experienced cracking problems due to excessive vibration. A solution is to 
Incorporate a well-designed friction dampers to attenuate blade vibration. The 
subject study, a cooperative effort between NASA Lewis and Carnegie-Mellon 
University, represents an application of recently developed friction damper 
technology to the SSME high-pressure oxygen turbopump. 

The major emphasis of study was the contractors design known as the two- 
piece damper. Damping occurs at the frictional Interface between the top half 
of the damper and the underside of the platforms of the adjacent blades. The 
lower half of the damper Is an air seal to retard airflow In the volume between 
blade necks. 

A bench test apparatus was developed to conduct an extensive set of exper- 
iments on the two-piece damper. The bench test apparatus was successful. The 
normal load was applied through a fishhook-pulley- weights system. The weights 
were varied from 0.25 to 70 lb. The blade support was tuned to simulate blade 
natural frequencies at pump operating temperatures and speeds (approximately 
9.5 KHz for the edgewise bending mode). The excitation system consisted of an 
electromagnet and a small chip of transformer iron mounted on the blade. Blade 
response, measured with a miniature accelerometer, was up to 700 G's tip accel- 
eration, 0.1-mll tip displacement, and 1 700- ps 1 stress at the crack location. 

Analytical models were fit to the experimental data and used to extrapo- 
late the results to pump operating conditions. These extrapolations show that 
the best thickness for the two-piece damper Is 0.047 In. and that the perform- 
ance can be Improved by reducing the width by 15 percent. Even with these 
Improvements, the performance Is poor. An example of predicted damper perform- 
ance (pump conditions) Is attached. The linear portion of the curve, emanating 
from the origin. Is where the damper Is locked and not functioning. The curved 
portion, beginning at an Input level of 3 lbf, is where the damper starts to 
slip. The optimum condition Is when the damper slips for approximately one- 
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half of the vibratory cycle. Since the stress required to crack the blade is 
estimated to be between 5000 and 10000 psi, we conclude that the damper is not 
working in the range of interest and is a poor design. 

Two other damper designs were considered. A design, known as the 
X-damper, was found to be a poorer performer than the two piece damper. A tip 
damper design was found to have excellent performance. 

Whirligig tests are inconclusive concerning two-piece damper performance. 
However, these results strongly imply that the parameter controlling blade 
vibrations is tip clearance. Large tip clearance results in high stresses. 

The hot-fire test results from MSFC are the only one which imply that the 
two-piece damper is effective in reducing stresses. Although this may be the 
case, the lack of cracked blades may be due to something other than the two 
piece dampers. For example, if the tip clearances on the blades in the hot 
fire test are relatively tight, the blades would have low stresses without 
dampers . 

In summary, the two piece and X-dampers are poor performers. The tip 
damper is a good performer and should be more actively pursued. Because of the 
unknowns associated with tip clearance, the hot-fire test should not be used 
as a justification to incorporate the two piece damper into production. Addi- 
tional testing should be conducted to understand the role of tip clearance on 
blade response. 
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DEVELOPMENT OF AN INTEGRATED BEM FOR HOT 
FLUID-STRUCTURE INTERACTION 

G.F. Dargush and P.K. Banerjee 
State University of New York at Buffalo 
Buffalo, New York 


The design of hot-section components such as those in SSME is 
generally accomplished by examining the solid and the fluid individually. 
Typically, the analysis that is required to support the design process is 
performed by separate groups using different techniques and computer codes. 
For a solid, a two-step approach is undertaken, a thermal analysis to 
determine the temperature profile using either a finite difference or a 
finite element method and the subsequent stress analysis by finite element 
method. 

The finite element method for stress analysis has reached a very high 
level of maturity and sophistication. Principal codes have comprehensive 
facilities for material and geometric nonlinearities, inhomogeneity, 
anistropy, as well as considerable flexibility, in the specification of 
boundary conditions. A similar development of boundary element is well on 
the way as a part of the host program and has already proven to be 
extremely attractive. and cost effective for both linear and nonlinear 
stress analyses problems. 

Computational fluid dynamics is, however, at a much less developed 
stage. The governing equations are highly nonlinear and coupled. Finite 
difference method dominates the scene primarily due to its simplicity and 
consequently one has to pay a price for this simplicity. Extremely dense 
grids are required near boundary and in regions of high velocity gradients. 
Varying requirements of discretization in the flow field, boundary layer, 
etc. can often lead to stability and accuracy problems. Finite element 
methods have also been applied recently, but it also suffers from the same 
restrictions. 

There is therefore a need to examine an alternative approach that 
considers the fluid flow and the solid deformation as a coupled problem 
using boundary element method. Since this method uses boundary equations. 
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it would be easier to develop a coupled analysis and since the 
nonlinearities in the fluid are primarily confined to regions, it may be 
possible to develop an efficient and cost effective solution. 

METHOD GP ANALYSIS 

Because the problem is essentially a coupled one, any separate 
treatment of each part (fluid and solid) would require the introduction of 
a complex set of artificially devised boundary conditions at the common 
interface to take care of the complex flow and heat transfer phenomenon. 
The problem is schematically depicted in Figure l, where the solid region 
is subjected to pressure (p), temperature (T) resulting in displacements 
(u) and stresses. The flow of fluid from both inside the solid and outside 
is governed by the velocity (v), pressure (p) and temperature (T) as well 
as the boundary displacement (u). 

A boundary integral representation for the ait ire coupled problem has 
been recently completed. The formulation is based on the fundamental 
analytical solution of the Navier-Stokes equation for the fluid velocity in 
an infinite domain. This fundamental solution was obtained by decomposing 
the Navier-Stokes equation into vorticity and dilation transport equations. 
A boundary integral formulation involving convolutions in time was then 
constructed in which the convective terms appear in the volume integral. 

The entire fluid formulation has been implemented and the linear part 
of. the viscous flow is already giving good results. It is hoped to get 
some nonlinear flew problems working within the next few months. Figure 2 
shows 4 boundary elements and l cell representation for the couette flew 
problem and Figure 3 shews that the exact analytical solution is reproduced 
by such a simple discretization. It is therefore likely that viscous flew 
at a high Reynolds number could be efficiently solved by this method. 
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FIGURE 1 

HOT FLUID STRUCTURE INTERACTION 


BOUNDRRY ELEMENT MESH 
COUETTE FLOW 



FIGURE 2 
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STERDY STATE COUETTE FLOH 

Test Problem - Incompressible Fluid 
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